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ABSTRACT 

The  purpose  of  this  report  is  to  present  the  results  of  the  first  phase 
of  a  detailed  analytical  study  of  the  combined  heat,  mass,  and  momentum 
transfer  processes  that  occur  between  a  liquid  stream  and  a  concentric, 
higher  velocity,  annular  gas  stream  when  both  the  liquid  stream  and  the 
annular  gas  stream  are  confined  in  an  axisymmetric  constant  pressure, 
converging  channel  or  an  axisymmetric  constant  area,  variable  pressure 
channel. 

The  study  began  with  an  effort  to  determine  whether  the  liquid  stream 
entering  the  r.  <cing  section  of  an  ejector-like  device  remained  intact  in  a  jet 
or  atomized.  The  conclusions  were  that  for  the  range  of  varibles  of  inter¬ 
est  the  probability  was  that  there  was  no  significant  jet  breakup.  Possible 
entrainment  of  liquid  by  the  concentric  gas  stream  has  been  neglected  at 
this  time. 

Having  establiiihed  an  analytical  model  based  on  the  assumption  that 
the  liquid  stream  remains  intact,  the  basic  equations  were  numerically 
solved  for  a  variety  of  conditions.  Using  the  axial  distance  required  to 
condense  the  vapor  from  the  concentric  gas  stream  on  to  the  liquid  jet  as 
a  measure  of  performance,  (the  shorter  the  better),  it  was  found  that: 

1.  A  constant  pressure  converging  mixing  section  is  more 
effective  than  a  constant  area  mixing  section. 

2.  Increasing  the  number  of  liquid  jets  (constant  total  jet  flow 
area)  shortened  the  mixing  section. 

3.  Increasing  the  Mach  number  cf  the  gas  stream  at  the  mixing 
section  inlet  shortened  the  mixing  section. 

4.  The  amount  of  noncondensable  present  in  the  gas  stream 
(from  20  to  50  percent  by  weight)  aid  not  significantly  affect  the  dis¬ 
tance  required  for  condensation  but  did  affect  the  shape  of  the  mixing 
section. 
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INTRODUCTION 


General: 

The  purpoee  of  this  report  is  to  present  the  analysis  and  results 
of  our  study  of  the  combined  heat,  mass,  and  momentum  transfer  pro> 
cesses  that  occur  between  a  liquid  jet  and  a  concentric,  higher  velocity, 
annular  gas  stream  when  both  the  jet  and  annular  gas  stream  are  con¬ 
fined  in  an  axisymmetric  channel.  This  particular  investigation  was 
undertaken  as  part  of  our  continuing  program  in  the  fundamental  study 
of  ejector  processes  for  the  Naval  Underwater  Ordnance  Station, 
Newport.  Rhode  Island. 

Generally  speaking,  most  ejector  type  devices  are  today  designed 
by  a  combination  of  empirical  data  and  one -dimensional  thermodynamic 
equilibrium  analyses.  The  results  of  such  an  approach  have  been  satis¬ 
factory  i;i  most  applications.  However,  when  two-phase,  two-component 
flow  is  involved,  e.  g. ,  a  liquid  water  jet  and  a  concentric  water  vapor- 
carbon  dioxide  gas  mixture,  this  traditional  approach  has  not  yielded 
devices  of  high  performance  nor  has  it  predicted  or  explained  all  of 
the  phenomena  actually  observed.  Apparently,  the  addition  of  a  non¬ 
condensable  gaseous  component  and  the  restating  diffusiqn  processes 
changes  significantly  the  heat  and  momentum  transfer  processes  in  the 
system. 

If  one  examines  the  schematic  diagram  in  Figure  1,  one  can  see 
that  the  ejector  devices  under  coniiideration  here  have  three  regions  of 
interest;  the  nozzles,  the  mixing  section  and  the  diffuser.  The  nozzle 
processes  are  well  understood  and  nozzle  performance  can  generally 
be  accurately  predicted  (exceptions  are  flashing  flow  in  the  liquid 
nozzle  and  condensing  flow  in  the  gas  mixture  nozzle). 

We  chose  to  begin  our  study  with  the  diffuser  section  of  the 
ejector  and  Reference  6  presents  our  initial  attempt  to  explain  some 
of  the  experimental  data  presented  in  Reference  3.  The  results  of  this 
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early  analytic^  study,  even  though  one-dimensional  and  without  consider¬ 
ation  of  rate  processes,  indicated  that  the  only  way  one  could  rationalize 
the  experimental  data  was  to  assume  that  the  two-|diase,  two-component 
flow  entering  the  diffuser  section  was  not  at  thermodynamic  equilibrium; 
i.  e. ,  the  gas  and  liquid  components  were  at  different  bulk  temperature 
levels  and  that  the  two  phases  were  flowing  at  different  average  velocities. 
The  next  step  was  to  try  to  predict  the  diffuser  section  inlet  conditions 
and  hence  our  present  study  of  what  could  be  the  mixing  section  process. 

Jet  Breakup  and  Atomization: 

Our  study  of  the  mixing  section  process  began  by  trying  to  estab- 
.  lish  an  analytical  model  of  the  processes  involved.  Of  major  importance 
was  the  question  of  whether  the  liquid  jet  entering  the  mixing  section 
remained  a  solid  jet  or  atomized.  Appendix  A  of  this  report  contains  the 
details  of  our  study  of  the  literature  on  atomization  and  jet  breakup.  The 
results  of  this  study  of  the  jet  breakup  and  atomization  literature  indicated 
that  there  was  a  good  probability  that  the  liquid  jet  did  not  break  up  to  any 
significant  extent.  (Neither  the  entrainment  of  the  liquid  nor  the  waviness 
of  the  liquid  jet  caused  by  the  scrubbing  action  of  the  gas  flow  has  been 
considered  at  this  time.)  In  addition,  visual  results  in  Reference  16  also 
report  no  significant  jet  breakup.  Thus,  as  a  starting  point,  we  have 
assumed  the  liquid  jet  does  not  break  up. 

This  conclusion  and  the  choice  of  our  first  analytical  model  have 
been  reached  on  the  basis  of  a  limited  amount  of  information,  for  in  spite 
of  the  large  amount  of  work  in  existence  on  atomization  (7-13)*  and 
turbulent  jets  (15,  16),  little  significant  quantitative  data  on  the  behavior 


*  Numbers  in  parentheses  refer  to  references. 
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of  a  high  velocity  liquid  jet  surrounded  by  a  higher  velocity  axisymmetric 
annular  gas  stream  was  found.  Therefore,  one  of  the  major  recommenda¬ 
tions  of  this  study  proposes  to  fill  the  gap  in  our  knowledge  of  these  funda¬ 
mental  processes  of  jet  breakup  and  atomization  by  further  study  in  this 
area. 

Mixing  Section  Analytical  Models: 

Along  with  the  derivations  of  the  conservation  equations,  Appendices 
B  through  £  of  this  report  present  the  details  of  the  computational  procedures 
for  the  three  analytical  models  considered  herein.  All  analytical  models 
considered  thus  far  assume  that  the  liquid  jet  remains  intact  with  ho  liquid 
entrainment  by  the  concentric  annular  gas  stream  and  no  liquid  jet  surface 
area  increase  due  to  waves  on  the  liquid  jet  surface.  Table  I  presents  a 
summary  of  the  three  mixing  section  analytical  models  and  indicates  the 
differences  between  the  models  while  Table  II  summarizes  the  assumptions 
common  to  each  of  the  mixing  section  analyses.  Table  III  is  a  summary  of 
the  input  data  and  range  of  variables  for  which  computations  were  carried 
out*.  In  addition  to  the  detailed  analyses,  derivations,  and  computational 
procedures  presented  in  Appendices  B  through  E,  certain  auxiliary  informa¬ 
tion,  such  as  data  sources,  correlation  sources,  etc.  ,  are  presented  in 
Appendices  F  and  G. 

Method  of  Solution; 

Generally  speaking,  the  method  of  solution  for  each  analytical 
model  is  to  simultaneously  solve  the  conservation  equations  of  mass, 
momentum,  and  energy  for  the  combined  heat,  mass,  and  momentum 


*  Because  of  space  limitations,  all  calculated  results  are  not 
presented  herein. 


transfer  processes  between  the  liquid  jet  and  the  concentric  annular  gas 
stream.  The  actual  computational  procedure  employed  in  each  analytical 
model  is  different  due  to  the  differences  in  the  models.  However,  the 
numerical  technique  employed  in  each  case  involves  a  trial  and  error 
iterative  procedure.  The  mixing  section  is  divided  into  a  number  of 
intervals  so  that  difference  type  equations  are  involved.  In  each  interval 
of  mixing  section  of  each  case,  the  solution  is  built  around  the  initial 
estimate  of  a  bulk  gas -liquid  jet  inte,rface  temperature  and  successively 
improving  the  initial  guess  until  all  rate  processes  and  conservation 
laws  for  each  interval  are  satisfied. 

The  iterative  technique  referred  to  above  is  an  extension  of  the 
method  presented  in  References  17,  18  and  19.  The  only  significant 
difference  between  the  work  referred  to  above  and  that  presented  herein 
is  that  changes  in  kinetic  energy  of  either  the  gas  or  liquid  stream  have 
been  included  as  significant  in  ciadculating  total  heat  transfers.  This 
modification  of  the  earlier  results  was  necessary  because  of  the  high 
relative  velocities  encountered  in  some  of  the  cases  considered  below. 

At  present,  the  analytical  models  herein  are  restricted  to  use 
of  input  data  for  which  the  condensable  vapor  state  at  the  mixing  section 
inlet  plane  is  approximately  a  saturated  vapor  state. 
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TABL£  I 

SUMMARY  OF  MIXING  SECTION  ANAJ^YTICAL  MODELS 

Analytical  Model  I: 

1.  constant  pressure  mixii^  ejection 

Z.  no  noncondensable  gas  in  gas  stream  (gas  stream  completely 
condensable) 

Analytical  Model  Ha: 

1.  constant  pressure  mixing  section 

2.  condensable  vapor  and  noncondensable  gas  present  in  gas 
stream 

Analytical  Model  Ilb: 

1.  diners  from  Analytical  Model  Ila  in  that  a  simplified 
relative  velocity  is  assumed  between  concentric  streams 

Analytical  Model  III: 

1.  constant  area  mixing  section 

2.  condensable  vapor  and  noncondensable  gas  present  in  gas 


stream 
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TABLE  n 

SUMMARY  OF  ASSUMPTIONS  COMMON  TO  ALL  MIXING 

SECTION  ANALYSES 


1.  The  flow  is  steady. 

2.  The  Bow  is  one -dimensional. 

3.  The  liquid  is  incompressible. 

4.  The  ejector  processes  are  adiabatic. 

5.  Height  changes  are  negligible. 

6.  Surface  tension  effects  are  negligible. 

7.  There  is  no  shaft  work. 

8.  The  liquid  jet  reinains  intact. 

9.  There  is  no  entrainment  of  liquid  ir  the  gas  stream. 

10.  The  liquid  and  gas  streams  are  concentric  with  the  gas  Bow  being 
annular. 

11.  The  vapor  enters  the  mixing  section  in  a  saturated  or  nearly 
saturated  state. 


SUMMARY  OF  REQUIRED  INPUT  DATA  AND  RANGE  OF 
INPUT  DATA  USED  IN  MIXING  SECTION  COMPUTATIONS 
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RESULTS 


Jet  Breakup  and  Atomization; 

The  calculated  results  shown  in  Figures  2  through  6  compare  the 
predictions  of  a  variety  of  jet  breakup  and  atomization  correlations  for  a 
range  of  possible  operating  conditions.  A  more  detailed  discussion  of 
these  results  and  the  correlations  from  which  they  were  obtained  is  given 
in  Appendix  A. 

Mixing  Section  Analyses: 

1.  Constaint  Pressure  Mixing  Section,  No  Noncondensable  (Model  I): 

Figures  7  through  10  show  some  of  the  calculated  results  for  one 
set  of  input  data  to  analytical  Model  I  plotted  versus  distance  from  the 
mixing  section  inlet  plane.  The  figures  are  generally  self-explanatory 
and  show  how  various  temperatures  (Figure  7),  various  velocities 
(Figure  8),  various  Reynolds  numbers  (Figure  9)f  various  heat  transfer 
coefficients  (Figure  9).  and  the  liquid  jet  diameter  and  mixing  section 
diameter  (Figure  10)  vary  with  distance  along  the  mixing  section.  The 
percent  of  total  inlet  vapor  flow  condensed  versus  the  distance  from  niix- 
ing  section  inlet  plane  is  shown  in  Figure  26. 

2.  Constant  Pressure  Mixing  Section,  with  Noncondensable  (Models  Ua  and 
lib): 

Figures  11  through  14  show  some  of  the  calculated  results  for  one 
set  of  input  data  to  Analytical  Model  Ila,  plotted  with  distance  from  mixing 
section  inlet  plane  as  the  abscissa.  The  figures  are  generally  self-explana¬ 
tory  and  show  how  various  temperatures  (Figure  11),  various  velocities 
(Figure  12),  various  Reynolds  numbers  (Figure  13),  various  heat  transfer 
coefficients  (Figure  13),  the  mass  transfer  coefficient  (Figure  13),  and  the 
liquid  jet  diameter  and  mixing  section  diameter  (Figure  14)  vary  with 
distance  along  the  mixing  section. 
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Figures  15  through  18  show  the  results  calculated  from  Analytical 
Model  Hh  for  the  same  set  of  input  data  used  in  Analytical  Model  Ha.  These 
figures  are  comparable  to  Figures  11  through  14,  and  show  die  differences 
resulting  from  the  assumption  of  a  simplified  relative  velocity  between 
phases  in  the  mixing  section. 

The  percent  of  total  inlet  vapor  flow  condensed  versus  the  distance 
from  mixing  section  inlet  plane  is  shown  in  Figure  26  for  Models  lla  and 
Ilb.  Also,  for  Model  Ila  the  percent  of  total  inlet  vepor  flow  condensed 
versus  the  distance  from  mixing  section  inlet  plane  is  shown  for  various 
inlet  Mach  numbers  (Figure  23),  for  various  inlet  flow  rate  ratios  of 
vapor  to  total  gas  mixture  (Figure  24),  and  for  various  niunbers  of  liquid 
jets  (Figure  25). 

3.  Constant  Area  Mixing,  with  Noncondensable  (Model  III): 

Figures  19  through  22  show  some  of  the  calculated  results  for  one 
set  of  input  data  to  Analytical  Model  III  plotted  with  distance  from  mixing 
section  inlet  plane  as  the  abscissa.  The  figures  are  generally  self-explana- 
toiy  and  show  how  various  temperatures  (Figure  19),  various  velocities 
(Figure  20),  various  Reynolds  numbers  (Figure  21),  and  the  liquid  jet 
diameter  and  mixing  section  diameter  (Figure  22)  vary  with  distance  along 
the  mixing  section. 

The  percent, iof  total  inlet  vapor  flow  condensed  versus  the  distance 
from  mixing  section  inlet  plane  is  shown  in  Figure  26.  Figure  22  shows 
how  the  total  static  pressure  in  the  mixing  section  varies  from  the  mixing 
section  inlet  plane. 

Comparis on  of  Mixing  Section  Analytical  Results: 

In  Figures  23  through  26  we  compare  how  the  total  percent  of  vapor 
condensed  (based  upon  amount  possible  to  condense)  varies  versus  the  axial 
distance  from  the  mixing  section  inlet  plane  for  a  variety  of  conditions. 
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Figure  23  shours  the  effect  of  increasing  the  gas  mixture  mixing  section 
inlet  Mach  number  while  Figure  24  shows  the  effect  of  increasing  flie 
amount  of  nonccmdehsable  in  the  inlet  gas  mixture.  Figure  25  shows  the 
effect  of  increasing  the  number  of  liquid  jets  (fixed  total  liquid  flow)  and 
Figure  26  shows  the  effect  of  changing  the  mixing  section  analytical 
model.  Insofar  as  possible,  in  each  of  Figures  23  through  26  all  the 
conditions  are  held  constant.  As  before,  each  of  the  figures  is  self- 
explanatory. 

The  most  effective  set  of  operating  conditions  (using  shortest 
distance  required  to  achieve  90  percent  condensation  as  a  criterion  of 
effectiveness)  is  Analytical  Model  Ilb  -  constant  pressure  mixing  section 
with  a  simplified  relative  velocity  between  phases.  Unfortunately,  this 
model  for  the  relative  velocity  between  phases  is  too  optimistic  -  the 
relative  velocity  model  used  in  Analytical  Model  Ua  is  more  realistic 
and  thus  Analytical  Model  Ua  with  eight  liquid  jets  is  the  most  effective. 

The  results  in  Figure  25  show  that  for  a  given  amount  condensed, 
the  distance  required  decreases  as  the  number  of  liquid  jets  increases. 

We  did  not  go  beyond  eight  liquid  jets  at  this  time  because  our  assumption 
of  no  significant  atomization  would  not  be  valid. 

Comparison  of  Figures  13  and  21  shows  why  a  constant  pressure 
mixing  section  is  more  effective  than  a  constant  area  mixing  section.  The 
Reynolds  numbers,  and  hence,  the  heat  and  mass  transfer  coefficients, 
have  a  higher  average  value  in  the  constant  pressure  mixing  section  than 
in  the  constant  area  mixing  section.  This  is  because  the  constant  pressure 
mixing  section  converges  as  condensation  takes  place  and  thus  the  average 
gas  mixture  velocity  remains  higher  than  in  a  constant  area  mixing  section. 
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CONCLUSIONS 

Jet  Breakup  and  Atomization: 

The  major  conclusiomi  of  the  jet  breakup  and  atomization  study  are 

that: 

1.  The  mechanics  of  atomization  of  a  liquid  jet  in  a  high 
velocity  gas  stream  are  not  well  understood. 

2.  The  available  correlations  indicate  that  in  the  range  of 
volumetric  and  mass  flow  rate  ratios  of  liquid  to  gas 
under  consideration  in  the  condensuctor,  the  predicted 
liquid  jet  breakup  lengths  are  large  compared  to  the 
condensuctor  mixing  section  length  and  the  final  mean 
drop  sizes  are  of  the  same  order  of  magnitude  as  the 
initial  liquid  jet  diameter  (secondary  stream  nozzle 
exit  diameter). 

3.  The  mixing  process  in  the  mixing  section  of  the  conden¬ 
suctor  most  probably  involves  heat,  mass,  and  momentum 
transfer  from  the  gas  phase  to  a  liquid  jet  with  no  significant 
breakup  or  atomization  of  the  liquid  jet  occurring. 

Mixing  Section: 

The  general  conclusions'  of  the  mixing  section  study  are  that: 

1.  A  constant  pressure  mixing  section  requires  a  shorter 
axial  distance  for  condensation  of  an  equal  amount  of  vapor 
than  does  a  constant  area  mixing  section  for  identical 
mixing  section  inlet  conditions. 

2.  Increasing  the  number  of  liquid  jets  with  constant  total 
liquid  flow  area  and  liquid  inlet  velocity  in  any  mixing 
section,  decreases  the  axial  distance  required  for  a 
given  amount  of  condensation. 


% 

i 


14 


3.  lacreasing  the  ratio  of  noncondensable  gas  to  condensable 
vapor  for  a  constant  total  gas  flow  in  the  constant  pressure 
mixiu^  section,  does  not  significantly  affect  the  distance 
required  to  condense  90  percent  of  the  vapor  flow  as  long 
as  the  mixing  section  has  proper  profile. 

4.  Increasing  the  gas  i^ase  mixing  section  iiilet  Mach  number 
in  the  constant  pressure  mixing  section,  decreases  the 
axial  distance  required  for  condensation. 
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RECOMMENDATIONS 

Jet  Breakup  and  Atomization: 

The  principle  recommendations  are  that  a  transparent  mixing  sec> 
tion  be  provided  for  future  testing  of  the  condensuctor  to  permit  visible 
records  of  the  mixing  process  to  be  obtained  by  means  of  high-speed 
photography.  This  will  make  it  possible  to  verify  available  correlations 
on  jet  breakup  or  develop  new  ones. 

Mixing  Section; 

It  is  recommended  that  the  mixing  section  analytical  study  be 
continued  in  order  to  relax  some  of  the  .restricting  assumptions  of  the 
present  analyses  and  to  include  the  diffusing,  or  diverging,  section  of 
an  ejector-like  device  as  a  continuation  of  the  mixing  section  process. 
Some  specific  recommendations  are: 

1.  Allow  for  liquid  entrainment  in  the  annular  gas  stream 
to  determine  if  two-phase  flow  compressibility  effects 
are  involved  and  to  determine  its  effect  on  axial  length 
required  for  condensation. 

2.  Allow  for  superheated  vapor  in  the  gas  stream  at  the 
mixing  section  inlet  in  order  to  extend  the  range  of 
variables  over  which  the  analysis  may  be  used. 

3.  Cut  off  the  mixing  section  at  some  point,  say  90  percent 
of  allowable  condensation  completed,  and  add  a  diffuser 
section  via  a  transition  section  in  order  to  have  a  com¬ 
plete  ejector  analysis. 

4.  Fix  a  mixing  section  and  diffuser  geometry  based  on 
some  desired  design  point  and  then  use  the  fixed  geometry 
to  examine  off- design  operation. 
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APPENDIX  A 

LIQUID  JET  ATOMIZATION  AND  BREAKUP 

Introduction: 

This  Appendix  is  concerned  with  an  investigation  of  the  possible 
mechanics  of  the  mixing  processes  occurring  in  the  mixing  section  and 
methods  of  predicting  the  degree  of  liquid  jet  breadcup  and  atomization. 

In  order  to  establish  an  analytical  model  for  a  Uieoretical  study  of  the 
combined  heat,  mass,  and  momentum  transfer  between  the  two  concentric 
streams  entering  the  mixing  section,  one  must  be  able  to  predict  the  length 
of  the  liquid  jet  at  breakup  and  the  mean  liquid  drop  sizes  resulting  from 
the  disintegration  of  the  liquid  jet. 

A  review  of  the  literature  on  atomization  and  liquid  jet  breakup 
was  made  to  determine  the  factors  which  affect  the  atomization  of  a  high 
velocity  liquid  jet  in  a  higher  velocity  concentric  gap  stream.  The  available 
correlations  for  liquid  jet  breakup  and  atomization  thus  found  were  used  to 
predict  liquid  jet  breakup  length  and  final  mean  drop  sizes  over  the  ramge 
of  variables  of  the  NUOS  condensuctor  experimental  test  runs  (Reference  3). 

Review  of  Literature; 

Beginning  with  Lord  Rayleigh’s  study  of  jet  instability  in  1868,  the 
phenomenon  of  liquid  jet  disintegration  has  been  the  subject  of  numerous 
theoretical  and  experimental  investigations.  Unfortunately,  most  of  these 
investigations  have  been  concerned  with  liquid  jets  discharging  into  a 
vacuum  or  into  ambient  atmosphere.  Almost  no  work  has  been  done  with 
liquid  jets  being  disintegrated  by  means  of  high  velocity  concentric  gas 
streams. 

Miesse  (11)  points  out  that  the  theoretical  and  experimental  analyses 
of  jet  disintegration  that  have  been  made  can  all  be  grouped  into  three 
general  classes,  or  into  a  combination  of  these  classes.  They  are: 
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"  1)  Small  Disturbance  Method.  An  arbitrary  small  disturbance  is 
imposed  upon  the  surface  of  a  liquid  jet  of  given  properties,  and  the  effect 
of  this  disturbance  is  studied  analytically  by  consideration  of  the  well 

established  physical  conservation  laws.  - Becaxise  the  phenomena 

considered  are,  by  nature,  nonlinear,  the  mathematical  solution  of  the 
problem  can,  in  most  cases,  be  effected  only  by  employing  the  linearizing 

assumptions  of  the  method  of  small  perturbations.  - An  obvious 

consequence  of  the  use  of  this  assumption  is  that  the  solution  obtained  is 
no  longer  (rigorously)  valid  when  the  magnitude  of  the  growing  disturbance 
exceeds  a  small  percentage  ---of  the  jet  radius. 

"2)  Method  of  Plaiisible  Conjectures.  For  the  phenomena  of 
secondary  atomization,  which  occur  subsequent  to  the  initial  growth  of 
the  disturbance,  mathematical  solutions  can  be  obtained  only  by  making 
plausible  conjectures  regarding  the  mechanism  of  the  governing  process. 
The  validity  of  these  conjectures  can  be  substantiated  only  by  the  degree 
to  which  they  effect  correlation  of  the  pertinent  experimental  data. 

"3)  Method  of  Dimensional  Analysis.  As  the  complex  combinations 
of  physical  processes  which  occuz*  in  the  disintegration  of  a  liquid  jet 
generally  defy  rigorous  mathematical  analysis,  the  pertinent  experimental 
data  can  frequently  be  correlated  effectively  by  consideration  of  the  dimen¬ 
sionless  groups  of  physical  varibles  which  have  an  effect  on  the  physical 
process. ” 

The  small  disturbance  method  has  been  used  to  study  jet  instability, 
wave  formation,  and  jet  breakup  by  R’ayleigh,  Tyler,  Weber,  Taylor,  and 
recently  by  Morrell  (12).  This  method  has  also  been  employed  by  Lane, 
Hinze,  and  Isshiki  (8)  to  predict  the  breakup  of  liquid  drops.  A  classic 
example  of  the  method  of  plausible  conjecture  is  that  of  Castleman  (13), 
who  utilized  Rayleigh's  theory  of  jet  breakup  to  postulate  ligament  forma¬ 
tion  as  an  intermediate  step  to  drop  formation.  Noting  that  liquid  droplets 
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seemed  to  be  torn  directly  from  the  main  body  of  the  liquid  jet,  he  showed 

that  ^^e  life  period  of  a  liquid  filament  from  which  the  drops  formed  is  in 

-5 

the  order  of  magnitude  of  10  seconds.  Although  Castleman  was  unable 
to  see  these  filaments,  recent  high  speed  photographic  techniques  have 
permitted  other  investigators  to  prove  their  existence.  The  method  of 
dimensional  analysis  tocorrelat  experimental  data  has  been  used  by  the 
majority  of  investigators. 

In  particular,  Ohnesorge  (as  described  by  Marshall  (7))  correlated 
the  experimental  data  of  Hanlein  by  plotting  the  Z-number  (ratio  of  Weber 
number  to  Reynolds  number  of  the  liquid  je^  versus  the  Reynolds  number 
of  the  liquid  jet  on  a  log«log  chart.  From  this  correlation  he  classified 
the  modes  of  jet  disintegration  in  four  groups  according  to  the  rapidity  of 
drop  formation  as  follows: 

"1)  Slow  dripping  from  an  orifice  without  jet  formation. 

”2)  Rayleigh  mechanism  of  jet  breakup  wherein  an  axially 
symmetric  disturbance  produces  breakup.  The  term  varicose  has  been 
applied  to  these  disturbances. 

”3)  Breakup  caused  by  disturbances  which  are  symmetrical  about 
a  helical  axis  starting  at  the  orifice,  as  treated  by  Hanlein  and  Weber. 

These  disturbances  are  sometimes  referred  to  as  sinuous. 

"4)  So-called  atomization  of  the  jet.  " 

"Since  for  a  given  liquid  and  orifice  size  the  Z-number  is  constant, 
a  variation  with  the  Reynolds  number  on  the  chart  follows  a  horizontal  line. 
Thus,  at  low  Reynolds  number.  Zone  I,  the  mode  of  breakup  follows  the 
Rayleigh  mechanism  and  as  the  Reynolds  number  increases  the  mode  passes 
into  Zone  II,  where  breakup  follows  a  lateral  motion  with  increasing  amplitude, 
and  the  jet  has  a  twisted  or  sinuous  appearance.  This'is  a  narrow  region,  and 
there  is  a  sharp  transition  from  this  zone  to  the  zone  of  atomization  at  the 
orifice,  Zone  III.  " 
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It  is  interesting  to  note  here  that  this  correlation  is  based  on  a  liquid 
jet  in  ambient  atmosphere.  If  the  liquid  jet  is  subjected  to  a  high  velocity 
gas  stream,  the  momentum  transfer  from  the  gas  phase  to  the  liquid  phase 
tends  to  accelerate  the  jet  as  it  progresses  through  the  gas.  Also,  acceler¬ 
ation  of  the  jet  tends  to  stabilize  it,  decreasing  the  amplitude  of  its  surface 
disturbances  and  keeping  it  intact. 

The  physical  properties  affecting  jet  disintegration  are: 

1.  liquid  nozzle  diameter,  D 

Lil 

2.  liquid  jet  velocity,  V 

XjI 

3.  relative  velocity  between  gas  and  liquid, 

4.  liquid  density,  p 

Xj 

5.  liquid  viscosity,  u 

Lj 

6.  surface  tension,  cr’ 


7. 

8. 
9. 


gas  density,  p 

8 

gas  viscosity,  p 

8 

geometry  of  mixing  chanber,  D 


Dimensionless  quantities  considered  in  various  analyses  are:. 

1.  Reynolds  number.  Re 

2.  Weber  number.  We 

3.  Z-number,  Z  =  We/Re 


4. 

5. 

6. 

7. 

8. 
9. 


density  ratio  of  gas  to  liquid,  p  /p^ 

8 

viscosity  ratio  of  gas  to  liquid,  p  /p.. 

8 

ratio  of  wavelength  of  disturbances  to  nozzle  diameter, 

ratio  of  breakup  length  to  nozzle  diameter,  L,  /D_  , 

b  Lil 

ratio  of  final  mean  drop  diameter  to  nozzle  diameter,  D  /D 

dr  Lil 

mass  flow  rate  ratio  of  liquid  to  gas,  oj 


10.  volumetric  flow  rate  ra.tio  of  liquid  to  gas,  w 


*  See  series  of  photographs  in  recent  EOS  Report  on  jet  condensers  (16). 
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Experimental  evidence  indicates  that  the  atomization  of  a  liquid  jet 
by  means  of  a  high  velocity  gas  stream  is  a  threefold  process  characterized 
by  the  following  distinct  phases: 

1.  The  liquid  jet  undergoes  surface  deformations  which  cause 
the  jet  to  break  up  when  the  amplitude  of  the  deformations 
become  sufficiently  large. 

2.  Liquid  filaments  are  detached  from  the  main  mass  of  the 
jet.  These  filaments  are  unstable  and  extremely  short¬ 
lived  (on  the  order  of  10"^  seconds)  and  immediately  form 
drops . 

3.  These  initial  drops  are  subjected  to  surface  deformation  and 
are  eventually  broken  down  into  successively  smaller  drops. 

Nukiyama  and  Tanasawa  (7)  and  Isshiki  (8)  have  investigated  liquid 
jet  breakup  and  liquid  droplet  deformation,  respectively.  Nukiyama  and 
Tanasawa  developed  equation  A9  to  predict  the  final  volume -surface  mean 
drop  size  for  the  overall  process  of  atomization.  Isshiki  derived  a  semi- 
empirical  expression,  equation  A13,  to  correlate  the  initial  and  final  drop 
size  for  the  third  phase  of  the  atomization  process. 

Both  of  these  investigators  include  liquid  to  gas  flow  rate  ratio  terms 
in  their  correlations  thereby  accounting  for  a  finite  gas  medium  and  the  trans¬ 
fer  of  momentum  from  the  gas  to  the  liquid  phase. 


Calculated  Parameters: 


Liquid  jet  Reynolds  number 
^Ll^Ll^L 


(Re)  = 


Liquid  jet  Weber  number 

^L^Ll 


a' 


1/2 


(Al) 


(We).  =  / 


/ 


(A2) 
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Liquid  jet  Z-r number 
Z  =  (We)./(Re)j 


(A3) 


Liquid  drop  Reynolds  number 

- ^7— 


(A4) 


Liquid  drop  Weber  number 

2  * 

P  V  ,  r 
/•ur-\  _  8 

(We)di  - 


(A5) 


Liquid  jet  breakup  length 


**(11)  Lj^  =  94.4D^j 


(7) 


Ll 

^L1 

cr' 

3,  /x. 


4  I  t  5 


"g 

1/2 


-I  1/8 


(12)  L, 


-z  . 

"'•I'’.  J 


1/2 


1/2 


(A6) 


(A7) 


(A8) 


*  Initial  drop  radii  [  (r  ,J  =  r_  , 
^  dr  max  LI 


**  numbers  in  parentheses  refer  to  references. 
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Final  mean  drop  radius 


(7)  r 


dF  2 


1410.  g* 


0.5 


+  191 


ra. 


0.45 


(1000  w  J 
vol 


1.5 


(A9) 


(7) 


(AlO) 


(10)  r^^  =  91.500. 


\  / 


IF 


rl 


(All) 


’^dF  "  ^L1 


gt 


1/3 


(8) 


-  \0.25 

^dFl 

^dl 


1.6 


23.5  +  0.000395 


"•^l'li^i 


(A12) 


0.5 


0. 25  2  P 

{We)di  l^L 


(Hw  )  C 


DI 


0.25 


0.  125 


In 


dl 


dF 


(We) 


dl 


(A13) 


Results; 

The  correlations  given  in  equations  A1-A13  were  used  to  calculate 
liquid  jet  breakup  length  and  final  mean  drop  radius  for  two  selected  experi¬ 
mental  test  runs.  Run  No.  1  represents  a  typical,  noncondensable  flow  case, 
i.  e.  ,  the  gas  phase  consists  of  a  mixture  of  water  vapor  and  noncondensable 
CO  gas;  Run  No.  2  represents  a  typical  condensable  flow  case  with  the  gas 

Li 

phase  consisting  entirely  of  water  vapor.  Test  data  for  the  two  cases  is 
listed  in  Table  IV  while  the  calculated  parameters  are  given  in  Table  V. 


fer;- 


A  computer  program  was  developed  to  perform  the  calculations  and 
is  set  up  to  perform  calculations  based  on  a  given  set  of  condensuctor  inlet 
conditions,  nozzle  exit  diameters,  as  well  as  the  static  pressure  at  mixing 
section  inlet.  In  order  to  investigate  their  effect  upon  the  calculated  para¬ 
meters,  three  of  these  conditions  may  be  varied  for  any  remaining  set  of 
inlet  conditions.  These  remaining  variables  are  mass  flow  rate  ratio, 
ratio  of  static  pressure  at  mixing  section  inlet  to  gas  stagnation  pressure 
at  ejector  inlet,  p^/p  and  ratio  of  liquid  stagnation  pressure  to  gas  Bta.g- 
nation  pressure  at  ejector  inlet,  This  computer  program  will  be 

useful  in  the  future  to  compare  calculated  values  of  liquid  jet  breakup  char¬ 
acteristics  with  experimental  test  results.  Such  a  comparison  can  establish 
the  accuracy  of  the  available  correlations  or  lead  to  more  accurate  ones. 

Figures  2  and  3  show  how  liquid  jet  breakup  length,  L^,  varies  with 
the  static  pressure  at  the  mixing  section  inlet  for  a  noncondensable  flow  case 
and  a  condensable  flow  case,  respectively.  The  actual  experimental  test 
pressure  is  indicated  in  each  figure.  The  correlations  given  by  equations 
A6  and  A8  show  that  the  length  of  the  liquid  jet  at  breakup  increases  as  the 
pressure  at  the  mixing  section  inlet  decreases.  This  is  caused  primarily 
by  the  decrease  in  gas  density  attending  the  drop  in  pressure.  The  resulting 
increase  in  relative  velocity  between  the  liquid  and  gas  phases  would  seem 
to  stabilize  the  liquid  jet  permitting  greater  lengths  before  breakup.  The 
correlation  given  by  equation  A7  does  not  account  for  the  effects  of  a  high 
velocity  gas  stream  surrounding  the  liquid  jet,  nor  the  initial  linu.fd  jet 
velocity.  Equation  A6  indicates  that  the  liquid  jet  will  be  longer  at  breakup 
if  the  initial  velocity  of  the  liquid  jet  increases . 

Figures  4  and  5  show  how  the  ratio  of  final  mean  drop  radius  to 

initial  jet  radius,  r,„/r-  , ,  varies  with  the  mass  flow  rate  ratio  for  a 

ux  Lil 

hone ondens able  fl,ow  case  and  a  condensable  flov  case,  respectively.  The 
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>  actual  test  value  of  mass  flow  rate  ratio  is  indicated  in  each  figure.  The 
results  of  four  correlations  given  by  equations  A9,  AiO,  AH  and  A12  are 
shown  in  each  figure.  Because  of  the  observed  transitional  nature  of  jet 
breakup  and  atomi.^ation,  it  is  highly  unlikely  that  atomization  occurs  at 
radius  ratios,  equal  to  or  greater  than  1.0.  The  range  of  re¬ 

sults  is  large  with  equation  A9  giving  a  fair  average  of  the  other  correla¬ 
tions.  Equations  A9  and  All  indicate  that  the  radius  ratio,  r  /rLl  , 
increases  as  the  mass  flow  rate  ratio,  u ,  increases.  This  trend  agrees 
^  with  the  results  obtained  from  the  correlations  on  liquid  jet  breakup  ‘ 

length  in  that  higher  initial  liquid  jet  velocity  tends  to  postpone  jet  break-  j 

up  and  droplet  formation.  Equations  AIO  and  A12  plot  as  straight  lines  in  | 

the  figures  with  values  from  equation  A 12  being  approximately  100  times 
greater  than  those  from  equation  AIO. 

Figure  6  shows  how  the  ratio  of  final  mean  drop  radius  to  initial 

jet  radius,  r,_/r_  , ,  calculated  from  equation  A9  varies  with  volumetric 
dr  Lil 

flow  rate  ratio  of  liquid  phase  to  gas  phase,  Three  points  are 

indicated  for  each  flow  case  at  various  values  of  static  pressure  at  the 

mixing  section  inlet,  p^.  This  plot  indicates  that  the  radius  ratio, 

r,_/r^  ,,  increases  as  the  volumetric  flow  rate  ratio,  w  ,,  increases. 
dF  Li  1  vol 

Equation  A13  may  be  used  to  calculate  final  mean  drop  radius  if 
the  initial  drop  radius  at  the  beginning  of  liquid  jet  breakup  is  known.  For 
purposes  of  these  calculations,  we  have  assumed  that  the  initial  drop  radius 
may  range  between  a  maocirnum  equal  to  the  initial  jet  radius  and  a  minimum 
based  upon  an  assumed  initial  drop  Weber  number,  equal  to  6.  5. 

Table  V  lists  the  initial  and  final  drop  radii  obtained  for  the  two  representa¬ 
tive  flow  cases.  In  both  cases  the  assumption  of  a  minimum  initial  radius 
does  not  produce  any  change,  whereas  the  assumption  of  a  maximum  initial 
radius  produces  a  growth  in  drop  size.  Thus,  the  correlation  of  equation  A13 
predicts  no  droplet  breakup  if  liquid  jet  breakup  is  assumed. 


t— iir** 


SUMMARY  OF  SELECTED  EXPERIMENTAL  DATA  USED 
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TABLE  V 

SUMMARY  OF  CALCULATED  RESULTS  FOR  JET  BREAKUP 
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APPENDIX  B 

CONST/ NT  PRESSURE  MIXING  WITH  NO  NONCONDENSABLE 

(ANALYTICAL  MODEL  I) 

Computational  Procedure: 

We  begin  by  calculating  the  following  quantities  at  the  ejector  inlet 
stagnation  state  (state  o)  and  at  the  entrance  to  the  mixing  section  (state  1). 


STATE  =  1. 

(Bl) 

T_  =T  /(T  /T^  ) 

Lo  go  go  Lo 

(B2) 

T  =  T 

LI  Lo 

(B3) 

Plo  =  Pgo  <PLo^Pgo' 

(B4) 

P  =  Pgo  (p/p 

go^ 

(B5) 

<b  ~  1 . 

^  o 

(B6) 

w  ,  =  w 

Vl  VO 

(B7) 

o 

II 

O 

(B8) 

(B9) 

w  =  w 
go  VO 

(BIO) 

w  .  =  w 

gl  go 

(Bll) 

(B12) 

V  -  c. 

8o 

(B13) 

^L1  ^dL 

1  / 

w  =  3600. 

Li 

Pl  '"li 

(B14) 

T  (R)*  =  T 
go 

+  459.69 
go 

(B15) 

T  (K)  =  5. 
go'^ 

T  (R)/9. 
go 

(B16) 

*  All  temperatures  are  in  degrees  F  except  as  indicated  by  (R)  or  (K)  which 
refer  to  the  Rankine  and  Kelvin  temperature  scales,  respectively. 
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CallSATT**  (p,  T  .  T  J 

go  sat 

T  =  T 
sat  1  sat 

Call  VAPTP  (T  ,  p  .  h  .  s  ,  v  ) 
go  ^go  go  go  go 

8  =  s 

gl  go 

Call  VAPTP  (T  p.  h  .  s  .  V  ) 

'  sat  ^  g  g  g 

If  (s  ,-s  )  B21.  B19,B26 
g^  g 

T  =  T  ^ 
gl  sat 

h  ,  =  h 
gl  g 

go  to  B27 


T  =  T 
gl  sat 

Call  SATLQ  (T^^.  v^.  s^.  h^) 

8.  =  8  -8- 
Ig  g  1 

''fg  ' 

h  ,  =  h .  +  X  h- 
gl  1  Ig 

go  to  B27 

CallSUPHTKT^^^j.  p.  s^. 
Call  VAPTP  (T^j.  p.  s^j. 


V  =  C ,  (2.  g  J  (h  -h  ,)) 
gl  dg  '  ®o  '  go  gl" 

h  ,  =  h  -Cj  ^  (h  -h  ,) 
gl  go  dg  'go  gl' 

If  (h  -h  )  B29.  B29,  B37 
gl  g 

T  =  T  ^ 
gl  sat 


1/2 


(B17) 

(B18) 


{B19) 

(B20) 

(B21) 

(B22) 

(B23) 

(B24) 

(B25) 

(B26) 

(B27) 

(B28) 

(B29) 


Fortran  Subroutines  are  identified  in  Appendix  G. 
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X  ,(R)  =  T  ,  +  459.69 
gl 

(B30) 

T^j(K)  =  5.  T^j(R)/9. 

(B31) 

X  = 

(B32) 

w  -  ST  X  w 

Vl  VO 

(B33) 

w  -  =  w  (1.  -X) 
cl  yo 

(B34) 

„  2  2 

w  (c  (T  -T  )  +  ^-  — ^  ) 
cl  y  pLi  gl  Ll  2.  g  J  / 

(B35) 

(B36) 

LI  "  Ll  C  -  w_ 

pL  L 

go  to  B40 

CallSUPHT2  p. 

(B37) 

T  ,(R)  =  T  ,  +■  459.69 
gl  gl 

(B38) 

T^^{K)  =  5.  T^j(R)/9. 

CallVAP(T^j,  T^l(K).  P.  (C/k)^.  W^,  a^.  J. 

(B39) 

k  C  ,) 

gl  Pgl 

(B40) 

pvl  pgl 

(B41) 

Pgl  Vgi) 

(B42) 

Dgi  =  (4.  (Agi+Aj^i)/n  '*■ 

{B43) 

CD  =  W  /W 

1  .  1^  gl 

(B44) 

CD  =  W  , 

o  1 

(B45) 

(w  /S.)  =  0. 

C  1 

(B46} 

1— ■ 

II 

o 

{B47) 

o 

11 

(B48) 
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• 

o 

II 

a 

(B49) 

phase  velocity  ratio  = 

(B50) 

bulk  temperature  ratio  * 

(B51) 

'  Sgi  “gi'Ni 

(B52) 

Call  VISLIQ  (T  . 

(B53) 

\l=  ^^gl-^Ll 

(B54) 

(B55) 

D^l=4.  A^l/(^(Dgi+c^D^l)) 

(E56) 

(Re)gi  =  3600.  D^,/(1Z.  p^,) 

(B57) 

T  =  T 
il  sat 

(B58) 

T.,(R)  =  T.,  +  459.69 

il  '  il 

{B59) 

T.j(K)  =  5,  T.^(R)/9. 

CallSATP  (T.^,  p^j) 

CallSATHFG  (T.^.  p^^. 

(B60) 

^fgc  1  ^gc  1  ^fc  1 

Call  SATLQV  (T.,,  v  ,) 
li  cl 

Call  VIS LIQ  (T.j, 

(B61) 

p  =  1.  /v  , 

^cl  cl 

{B62) 

V  =  0.  75  V  , 
gl  rl 

(B63) 

VL  =  0.25  V  , 

LI  rl 

(B64) 

E  =  0. 

(B65) 

(Re)^j  =  3600.  PlV'liDl,/(12. 


(B66) 
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(Re)^^j  =  3600.p^jV.^jD^,/(12. 

If  ((Re)_  -  -  12,500.)  B68,  B68,  B69 
X^l 

h,  .  =  1.86  (12.  k,  /D,  ,)  (Re),  (Pr).  (f‘n/(lcl>'’' 


(B67) 


LI 
go  to  B70 


L'  LI' 


'Ll 


'Ll 


"li  =  0  023  (12.  kj^/D^j)  (Re)j^j®  «  (Pr)^j®-'‘ 

Call  CASP4  (T..(K).  0^.  (€/k)  .  W  ,  pt 

ll  V  V  V  '^gll 

IF  ((Re)^^j  -  9000.)  B71,  B71.  B72 

»G,  l-*‘  (‘2-  >'g,/°el>  <^*>rgl‘^'  <»'gl/»‘gil)“' 

go  to  B73 

hci  =  0  ''gl'°el>  <^*>rgl“  ®  <*“^>gl""  ‘^l/Mgu)"' 

If((Re)rgi  -  1000,  )  B74,  B74,  B75 
f.  =  16.  /  (Re)  , 
go  to  B76 

f.  =  0.0014+  (0.  125/(Re)  °*  ^^) 

1  rgl 

T,=  3600.^f.p^j  V'^j^/(2.  g') 

<-JL'®i>l  =  ^1  f^il-'^Ll' 


(B68) 

(B69) 

(B70) 


(w  /S.)  = 

c  1  1 


'0L/®i>i  -  ''Gi<\r'^ii' 


2  2 
V  -V 

h.  ,  +  C  _  (T.,-T_  .)  +  -g.L..  ■  hlr. 
fgcl  pL  il  LI  2.  g^J 


T  =  3600.  (w  /S.),  V  ,/g' 
m  c  1  1  gl  ® 


(B71) 

(B72) 

{B73) 

(B74} 

(B75} 

(B76) 

(B77) 

(B78) 


(B79) 


T  =  T.  +  T 
1  m 


(B80) 
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LI  “ 

Pi. 

1  ‘^*>rglj 

=  V 

V,  ,  -  V'  , 

gl 

gl  ’ 

LI  LI 

^  ^  r  /t>  \  ^  I 

2.  1  8 

V  + - - 

8^  (0.046)  3600.  ^p 


1/2 


If  (E-50.)  B83,  B83.  B86 


(B81) 

(B82) 


“( 


<”c/Si)l  -  ("c/Sj) 

(Wc/Sj), 


(w^/S.)  =  (w^/Sj), 

E  =  E  +  1. 
go  to  B66 


B86.  B84,  B84 


(B83) 

(B84) 

{B85) 


h.  ,  +  C  _(T.,-T_  ,)  + 
fgcl  pL  il  LI 


2  2 
\l 


AT 


1 


Ui  =  (q^/S.)i/AT, 


(B86) 

(B87) 

(B88) 


At  this  point  we  have  calculated  all  variables  at  state  1.  We  now 


commence  our  calculations  for  the  first  state  2. 


Ca'O. 


STATE  =  STATE  +  1. 


AT  =c,AT 
g  1  g 

T  ,  =  T  ,  -AT 
g2  gl  g 

If  (T  -  -  T  J  B93.  B94,  B94 
gZ  sat 

T  =  T 
g2  sat 


(B89) 

(B90) 

(B91) 

(B92) 

(B93) 


I 
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<'*c>1e“‘=4''vo 

T  _(R)  =  T  ,  +  459. 69 
g2  g2 

T^2(K)  =  5.  T^^{R)/9. 

If  (0.99  -  c^)  B98.  B98,  BlOl 

If  (1.399  -  c^)  BIOO.  BlOO,  B99 

(w  )-_  =  0.  02  w 

c  12  VO 

go  to  BlOl 

(w  ),_  =  0.  01  w 

c  12  VO 

w  _  =  w  -  -  (w  )  - 
v2  vl  c  12 

w  -  =  w  ,  -  (w 
g2  gl  c'12 

T  =  T 
12  sat 

Ti2(R)  =T,2  +  459.69 
T.2(K)  =  5.  T.^  (R)/9. 

CallSATP  (T.^, 

CallSATHFG  (T.^,  P^2'  ^fc2> 

^fgc2  ’  ^gc2  ■  ^fc2 

^c2=^cl+<Vl2 

CallVAP(Tg2'  ^v'  %*  V'  V’ 

k  C  ,) 
g2  Pg2 


C  =  C  , 
pv2  pg2 

^^^^g2  *  Sg2  ^g2^^g2 


(B94) 

(B95) 

(B96) 

(B97) 

(B98) 

(B99) 

(BIOO) 

(BlOl) 

(B102) 

(B103) 

(B104) 

(BIOS) 


(B106) 

(B107) 


(BIOS) 


(Bi09) 
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Ah  =  C  ,  T  ,(R)  -  C  ,  T  _(R) 

V  pvl  gl  pv2  g2' 

(q  )i,  =  w  -  Ah 

V  12  v2  V 

V  =  V 
L2  LI 

V  =  V 

g2  gl 

V  -  X  V  - 

g2  gl 

„  „  ,  *  Viz  ^VlZ  *Ngcl  “^pvZ  •’’’gz'^il’* 

-I  ^  g_  \ 


L2  LI 


<ViZp  =  “- 


SL<*L  +  *el> 


E  =  0. 


(BllO) 

{Bill) 

(B112) 

(B113) 

(B114) 

(BUS) 


(B116) 


(B117) 

(B118) 


(w  /S.)  =  (w  /S.)  +  2. 

Call  SATLQV  (T.^, 

Call  VISLIQ  (T.^. 

"cZ  '  '''cZ 

ATgc  =  -  Tj,  -  T.^)/Z 

aT"  =  (C  ,  +  C  -)  at  /2. 

VC  pvl  pv2  gc 
^fgc  "  ^^fgcl 

^cL  ""  SL^^^il‘*'’^i2?^^*  ■  ^L2^ 

(q  ),,  =  w  ,  (Ah  +  AV  ) 

^v  12  v2  V  g 

(q  ),,  =  (w  ),-  (h.  +Ah  +Ah,+AV_) 
c  12  c  12  fgc  VC  cL  gL 


(B119) 

{B120) 


(B121} 

{B122) 

(B123) 

(B124) 

(B125) 

(B126) 

(B127) 

(B128) 

(B129) 


{B130) 
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(q^)12  =  (w^  +  w^j)Ayj^ 

^12  “  ^^^2  ^^€^2  ■‘'  ^*^1.^2 

Q, 


T,  , .  =  T,  ,  +  — 
L2  LI  C 


12 


pL  **li  '*'  *cl* 


=  V  -  ""lZ 


M 


In 


^^gZ-’^LZ' 


Call  VISLIQ  (Tj^. 

‘*■'’1  2  '  Sl  '‘lZ^’^L 


(w^/S.)  = 


(w^/S.),  -  (w^/S.)^ 

^  V  KTsT),  j 


c'  i'2 


<®ic>lZ  =  ‘"c’lz'  <’'c/®i> 

Sense  Light  1 

A  s  144.  w  /(3600.  p  V  ) 
g2  g2  '^gZ  g2' 


w. 


w 


“^L2  ” 


c2 


^L^L2 


^c2^L2 


^L2*<^*  ^L2'  ^2» 


1/2 


1/2 


/ 


3600. 


°g2  =  <^'  <^2"^L2>/'> 

=  4.  A^^/(ir  (D^2  +  -2  ^L2» 


^2" 


V  -  V 
g2  L2 


(B131) 

(B132) 

(B133} 

(B134) 

(B135) 


(B136) 

(B137) 

(Bi38) 

(B139) 

(B140) 

{B141) 

(B142) 

(B143) 

(B144) 

(B145) 

(B146) 
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(B147) 

(B148) 

(B149) 

(BiSO) 

(Bi51) 

(B152) 

-  (B153) 

(B154) 

(B155) 

(Bi56) 

(B157) 

(B158) 

{B159) 

(Bl6p) 


(3161) 
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L2 


■  ( 


(Re) 


L2 


(Re) 


rg2 


.  0.2 

2.  T  (Re)  g 

V  + 

(0;046)  3600.^p 


] 


1/2 


V’  =  V  ,-V',  , 

g2  g2  L2  L2 

If  (Sense  Light  1)  B139.  B164 


(B162) 

{B163) 


Call  GASP4  (T._(K).  O  ,  (^/k)  , 

i2  V  V  V  gi2 

If  ((Re)  ^-9000.)  B165.  B165.  Bi66 

rg2 

=  1.8MU.  (Re),g2*'^Pr) 


(B164) 


go  to  Bl66a 


0.8 


hpj.  0.027  (IE.  k^,/D^^) 

If  , ((Re),  --42,  5.00;.:)  BI67,  B167.  B168 


1/3 


go  to  B16.9 


0.8 


hj^2  *  0-  023  (12.  1Cj^/Dl2>  ^^®^L2'^’ 
«*l/®1>2  =  ^L2 


(w  IS.), 

c  1  2 


'  ^G2  ^-^g2~^i2^ 

V  ^-V 

K  ,+C  _  (T.-  -T^  ,)  +  ^ 

fgc2  pL  i2  L2  2.  g^J 


1/3,  ,  ,0.14 

'gZ  •'*g2^*‘gi2' 

(B165) 

,1/3,  ,  ,0.14 

■’gz  ‘V^^giz’ 

(B166) 

(B166a) 

1/3,  ,  ,0.14 

.Z  <''l2'''c2> 

(B167) 

0.4 

L2 

{B168) 

(B169) 

~2 

(B170) 

(q  /S.),,  =h_._(T  _-T.,V+  (w  /S.)_  (l 
’g  i  2  02  g2  i'2;  c  i  2  \ 


h.  ,+C  JT._-T_-)+  ,  _ 

fgc2  pL  i2  L2  2.g  J 


V  -V  N 
L2  \ 


(B171) 

(B172) 
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nwsT — 


In 


V(w^/Si)J 


h2  = 


iz.(Si)i,  ^ 

2 

^°L2‘^L1  ' 

2  " 

iw  °L2'^°L1 

^'=2  2.  ^ 

i  ‘  J 

_ 

“1  1/2 


K  (E-50.)  B177.  B177.  B179 


If 


<®i>l2-<Si>12o 

‘®i'l2 


-  ^2  I  B178.  B178 


go  to  B120 


^2  °  ^1  ■*'  ^12 
“2'<'"L''*c2*''"g2 

♦2=*v2/V 

phase  velocity  ratio  =  V  ,/V,  , 

g2  L2 

bulk  temperature  ratio  =  T  ,/T,  . 

g2  L2 

®i2=®ir+(Si'l2 


Q  =  Q  +  Q 
^2  '^1  '^12 


We  have  now  finished  calculations  for  state  2.  If  we  satisfy  the 
following  tests,  the  calculation  is  complete;  if  not,  we  reset  several 
variables  to  establish  state  1  as  the  beginning  of  the  next  mixing  section 
interval  and  recalculate  state  2, 


(B173) 

(B174) 

(B175) 

(B176) 


(B177) 

(B178) 

(B179) 

(B180) 

(B181) 

(B182) 

(B183) 

(B184) 

(B185) 
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If  (T  _-T_  -)  B212.  B212.  B186 

g4  L,C 


If  (w  -w  _)  B2I2.  B212.  BI87 

VO  c2 

(B186) 

If(1.399-c^)  B212.  B212.  B188 

(B187) 

^fgcl  ~  *‘fgc2 

(B188) 

pvl  pv2 

(B189) 

Sgl  ”  ^pg2 

(B190) 

^1=^2 

(B191) 

(Wc/S.)i=(Wc/S.)2 

(B192) 

^vl  “  ^v2 

(B193) 

T  =  T 

LI  L2 

(B194) 

T  s  T 
gl  g2 

(B195) 

Tgi(R)  =  T^2(R) 

(B196) 

V  =  V 
gl  g2 

{B197) 

V  =  V 

LI  L2 

{B198) 

V«  3  V' 

LI  L2 

(B199) 

V  =  V 
gl  g2 

(B200) 

T  =  T 
il  i2 

(B201) 

w  =  w  , 
cl  c2 

(B202) 

^gi  =  V 

(B203) 

^L1  ■  °L2 

(B204) 

D  =  D 
gl  g2 

(B205) 

L  =  L 

1  2 

(B206) 

(B207) 
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«1 

(B208) 

(Re)^l  ,  (R.)^3 

(B209) 

(B210) 

"gj "  V 

(B211) 

go  to  B90 

read  input  data  for  new  case. 


(BZ12) 
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APPENDIX  C 

CONSTANT  PRESSURE  MIXING  WITH  NONCONDENSABLE 
(ANALYTICAL  MODEI^  Oa  AND  Ilb) 

Computational  Procedura  -  AmJ-ytical  Model  Ha; 

We  begin  by  calculating  the  following  iquantities  at  the  ejector  inlet 
stagnation  state  (state  o)  and  at  the  entrance  to  the  mixing  section  (state  1). 


STATE  =  1. 

(ei) 

T^  =T  /(T  /T^  ) 

Lo  go  go  Lo 

(C2) 

^L1  *  ^Lo 

(C3) 

^Lo  ”  ^go  ^^Lo^^gb^ 

(C4) 

P*P  (p/p  ) 

go  go 

(C5) 

a  s  w  /  (w  +w  ) 

^  o.  VO  VO  nc 

(C6) 

w  ,  =  w 

Vl  VO 

(C7) 

♦  1  '*0 

(C8) 

w  !S  W  +  W 

go  nc  VO 

(e9) 

w  ,  ss  w 
gl  go 

(CIO) 

(Cll) 

=  ^dL  «o 

(C12) 

w^  =  3600.  V^j  Pj^Aj^j/144, 

(CIS) 

T  (R)#  =  T  +  459.69 
go  go 

(C14) 

T  (K)  =  5.  T  (R)/9. 
go  go 

(C15) 

^  All  temperatures  are  in  degrees  F  except  as  indicated  by  (R)  or  (K) 
which  refer  to  the  Rankine  and  Kelvin  temperature  scales,  respectively. 
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Call  SATT*<(p.  T  ,  T  J 

go  sat 

K  =1.3 
g 

Call  GASP  (T  (K).  w  .  w  .  W  .  W  ,  J.  C  .  K  ) 
go  nc  VO  nc  v  Pgo  go 

K  -1. 


(C16) 


gl  go  go 
T^j(K)  =5.  T^j(R)/9. 

Call  GASP  (T  ,(K).  w  .  w  W  .  W  .  J.  C  K  .) 

gl  nc  vl  nc  V  pgl  gl 

=  K 
go  g 

=  (K  +  K  ,)/2. 
g  go  gl 


If  (  K  -K  -  C  .)  C21,  C17.  C17 
I  g  go  I  4' 


V,=C,  Q.g  J(C  T  (R)  -  C  ,T 
gl  dg  L  o  pgo  go'  pgl 


o  pgo  go 


,j(R))] 


=  LSgo  '^go'^>  -  %-  “^pgo-fgo'W  -  Sgl^gl'^'^l'^pgl 

T  (K)  =  5.  T^j(Rj'/9. 
w  ,/W 

X  -  vl  V 

vl  ”  w  ,  w 

vl  nc 

w  W 

V  nc 

Call  CASPl  (T  ,(K).  O  ,  a  ,  a,  t*,  (C/k)  ,  (C/k)  W  ,  W  , 
gl  nc  V  V  nc  nc  v 

p.  X  ,,  ,  D'  ,  u  a  ,,  U  J 

^  vl  gl  1  ^vl  ^ncl 


(C17) 

(C18) 


(C19) 

(C20) 


(C21) 

(C22) 

(C23) 

{C24) 


(C25) 


Call  CASK  (p.  X^j,  J,  Tg,(R).  T^j  (K), 

S'  ,  S'  ,  cr  p  T  ,(K).  C  e  c  ,) 
nc  V  gl  '^gl  gl  pgl  pncl  pvl 

R  =  1545.45/({1. -X  ,)  W  +  X  ,W  ) 
gl  vl  nc  vl  V 


(C26) 


**  Forti*an  Subroutines  are  identified  in  Appendix  G. 
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=  <*0^1  Tgl 

(C27) 

(C28) 

Agi  =  1«.  Wgi/(3600.  Pgj  V^j) 

(C29) 

’’gl  '  <^1  ^ 

(C30) 

W  .  =  W_  /w  - 
1  L  gl 

(C31) 

u  ~  u  , 
o  1 

(C32) 

(w  /S.)  s  0. 

C  1 

(C33) 

*01  '  “• 

(C34) 

h  =  o. 

(C35) 

S.,  =  0. 
ll 

(C36) 

Qj  =  0. 

(C37) 

phase  velocity  ratio  *  ^gX^^Ll 

(C38) 

'^■''gi  =  '^Pgi  -^gi^i 

(C39) 

Call  VISLIQ(T^j.  #1^^) 

^1=1  V-'^LlI 

(C40) 

V  ,  =  0.75  V  , 
gl  rl 

(C41) 

V'i^,  =  0.25 

(C42) 

(C43) 

D^l=4.  A^l/WD^j+c^D^,)) 

(C44) 

{Re)^j=  3600.  p^jVgjD^,/(12.p^j) 

(C45) 

'®='gl  =  V^<Pgl°'l’ 

(C46) 

(C47) 

E  =:  0, 


(C48) 
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">• 

Aqi3  =  0. 

T.,(R)  =T.,  +  459.69 

ll 

T.j(K)  =  5.  T.j(R)/9. 

CallSATP(T.j. 

^vil  "Pel 

Call  SATHFG  (T...  p  .  h  h.  ,) 
ll  cl  gel  fcl' 

^fgcl  ^gcl  ^fcl 

Call  VISLIQ  (T. , .  u  J 
ll  ^cl 

I£({Re)^j-12,500,)  C56,  CS6,  C57 

“li  '  ‘  (12-  V°Ll' 

go  to  C58 

Hl,  =  »•  (12-  (P-in®'" 


(Rfe) 


^^1=  3600.  D^^/(12. 


CallGASP3  (T.,(K),  <7  ,  a  ,  a,  t*.  (c/k)  .  (e/k)  .  W  , 

ll  nc  V  V  nc  nc 

W  ,  p,  p  U  .,) 

V  vil  '^gil 

If  ((Re)  ,  -  9000.)  C59,  C59,  C60 

rgl 

"gI  =  (*'•  <P')gl‘''  <'‘gl''^gll>‘’' 

go  to  C61 

hoi  =  0-027  (12.  (R.)^^/-0  (Pr)^//'  (Mgi/dgu*"' 


(K^'P'  ),  = 

G  nc  ’ 


i=^Gi  (^gi'^'iV 


1/3 


(C49) 

(C50) 

(C51) 

(C52) 

(C53) 

(C54) 

(C55) 


(C56) 

(C57) 

(C58) 


(C59) 


(C60) 


(C61) 
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Call  SATLQV  (T.,,  v  J 
jlJL  cX 


Pc,  =  1- 

1 

cl 

T  =  3600.  (w^ASO  V  /g- 
m  c  1  1  gl 

V  = 

LI 

Pgl.  /  («*>L1  f 

[Pl.  \<^')rgl/ 

V  .  = 

gl 

1  I 

P*  s:  ■ 

^vl  ^  ^vil 

ncl 

14.  696  In  ( 

K  =(K  -P'  )  /P« 

Gl  G  pc  1  ncl 

2.  T  (Re).  ®  V 

V'  ^ 

ffl  2 


(0.046)  3600.  p 


L  J 


1/2 


(Wc/Si)i  (X^j  P-P.,,V14.696 


vir 


‘g'  i'l  or  gi  ii 


fgcl  '  pL'*il  "LI' 


V  ^-V 

tV  LI 

2.  g^J 


(C62) 

(C63) 

(C64) 

(065) 

(C66) 

(C67) 

(C68) 

2 


(C69) 

(C70) 

Aqj  =  (q^/Sj),  -  (qg/S.)j 

{C71) 

^ill  '  ^U2 

(C72) 

Aqil=Aqi2 

(C73) 

If  (E-50.)  C74,  C74,  C82 

I£(  !  (T.,-T.j2)/T.j  I  -  Cj)  C82,  C75.  C75 

(C74) 

T  =  T 

il2  il 

(C75) 

Aqi2  =  Aqj 

(C76) 

E  s  E  +  1 . 


(C77) 
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If  (E-3.)  C78.  C79,  C79 


go  to  C80 

^il  “  ^^il2  ^^ir^ill  ^^12^  ^ 

(G79) 

If  (T  ,-T.,)  G81V  €81,  C51 
=  gl  il 

(G80) 

T.,s:T  ,-'0. 1 

il:  gl  1 

(G81) 

go  to  C51 

1  gl  LI 

(C82) 

(G83) 

.  ^mperafure  ratio  s 

(C84) 

At  this  pbint  we  have  calculated  all  variables  at  state  1.  We  now 

^commence  pur  calculations  for  the  first  state  2» 

STATE  r  SI  ATE  +  d. 

(C85) 

AT  =c,.AT 

8  ^  ? 

(C8b) 

If  (AT  -  C.-).  G88,  C87,  G87 
g  <: 

AT  =  f_ 
g  2 

(C87) 

T  ,  =  T  ,-AT 
g2  gl  g; 

.(G88) 

If  (T  ^-t  J  G90.  G89,  C89 
g2  sat 

T  .  =  T  ^  -5. 
g2  sat 

(C89) 

if(Tg2”'^LrV 

(C90) 

%Z  °  ■•■  *6 

(C91) 

Call  SATP  (T^j, 

(C92) 

^g2  ”  ^  ^v2 

(C93) 
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^  ^  p.  V  /p 

-.'y*  -yi- 

C9h  C.9i.  C92 

T  =  T  -  1 
g2  g2  "• 

go  to  C92 


(C94) 


(C95) 


w 


v2 


W  W 
nc  V 


W 


nc 


w  _  =  w  +  w  , 
g2  nc  v2 

<Vl2' W 


w  =  w  -  W  _ 
c2  VO  v2 


u(M 

\W 

VO 


-  €3]  ClOO,  G234,  C234 


='^g2''‘  '*59.69 
Tg2(K)  =  5.  T^2(R)/9. 

Call  GASPl  (T^2<K).  or.  t*.  (c/k)^,  (f/k)^^,  W^. 

P’  \2’  %2’  ^2’  ^v2’  ^nc2’  ^*g2^ 
eallGASP2  (p.  W^.  J.  T^^(R). 

^  nc'  ^’v'  ^g2'  Pg2'  ^g2<^>'  Sg2*  Snc2'  Sv2^ 

<P'>g2  =  Sg2 

Ah  =  C  -  T  ,(R)  -  C  -  T  ,(R) 

V  pvl  gl'  '  pv2  g2'  ' 

(q  )i,  =  w  -  Ah 
V  12  v2  V 

Ah  =  C  -  T  ,  (R)  -  C  ,  T  ,(R) 
nc  pncl  gl  '  '  pnc2  g2^  ^ 

(q  =  w  Ah 

nc  12  nc  nc; 

'®"’g2  =  ®'2> 


(C96) 

(C97) 

(C98) 

(C99) 


(ClOO) 

(ClOl) 


(C102) 

(Ci03) 

(C104) 

{C105) 

(C106) 


(C107) 


(C108) 


■V,  =•  V.. 


V  -  «  V  - 

g2  gl 


T  =  T  + 
■  L2  lil 


(q  )i,  +  (q  +  (w  (h,  .+C  ,  (T 

nc  12  V  12  c  12  fgcl  pv2  g2  il 


E  =  0. 


^iZZ  =  “• 

^<>22  '  «• 

v,  *  v  , 

g2  gl 

^'l2  ^ 

<*c'®i'2  '  -  ^• 

t.  j(R),'=  T.^  +  459. 69 

T.j,(K)  =  5.  T,2(R)/9. 

Call  SATP  (T.^,  p^.^) 


*’c2  '  ’’viz 

Call  SA9’«rG  (T.^, 

CallSATLQV  IT.,,  v  ,) 
i2  c2 

'^.all  VJtSLia  (T._,  4 

i2 


Pc2  = 


c2 


^fgc2  ”  ^gc2  ^fc2 

AT  =  (T  ,  +  T  ,  -  T.,  -  T.,)/2. 
gc  gl  g2  il  i2' 

Ah  =  (C  ,  +  C  ,)  AT  II. 

VC  pvl  pv2  gc 


fgc  '  fg 


(C109) 

(Cl  10) 

(Clll) 
(Cl  12) 
(C113) 
(C114) 
(C115) 
(Cl  16) 
(C117) 
(C118) 
(C119 

(C120) 


(C121) 

(C122) 

(C123) 

(C124) 


(Cl  2  5) 
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i- 


'  '^Ll'-'^Lz'*  !  <^-  *oJ> 
^\l.  -  ('^gl'-^Lz'^  '  '^-  *o"' 


Ah  =C  fliilZi?  -T 
cL  pL  \  2.  L2 

(q  ) , ,  =  w  (Ah  +  AV  ) 

V  12  v2  V  g 

(q  )  1 ,  =  w  (Ah  +  AV  ) 

nc  12  " 


) 


nc 


nc 


g 


“■c'lZ  =  ('"'chz  *  \c  ^  ^•'cL,-^  ^'^gL> 

<'’l>12  '  <*L  +  *cl> 

*  ***nc’l2  ■'■  *%*12  ■'■  *‘*ch2  ■'■  *‘*L*12 

Q, 


T  »  T  +  — 
L2  ^L1  C 


‘12 


pL  <*L  ^  *cl> 


AT  =  T  -T 
2  g2  L2 

■  (T  -T  )  -  (T  -T  ) 

at  ^ ^ 

M  /  T  .-T, 


In 


(  gl  Ll\ 

^  ^g2-^L2  I 


Call  VISLIQ 

<*'’^>L2  *  Sl  ''l2/'^L 


(w^/S.)  = 


,  /‘*c/®i>l\ 


<®ic>12  =  <'*'c'l2/<*c/®i> 

Sense  Light  1 

\z  *  '"g2^'^^“°-  Pg2V’ 


(C126) 

(C127) 

(C128) 

(Ci29) 

(C130) 

(C131) 

(C132) 

(C133) 

(C134) 

(C135) 

(C136) 

(C137) 


(C138) 

(C139) 

(C140) 


4 

4' 


: 


(C141) 


w. 


w 


c2 


PtV 


L  L2 


^c2^L,2 


3600. 


°L2  =  ''l2/<^"2>) 


1/2 


1/2. 


-’  <°gl  ^  ^g2>  ^  ^• 


°e2*"-  ''l2» 

^r2  '  I  V'^LZ  I 


h2  = 


f. 

‘^•<Vl2  ' 

2 

'®L2-^Ll\  ' 

D_  _  +  D_  . 

'■  ( 

. 

/ 

Vg-(v^,.v^,)/z. 

V'g  =  (Vg,  ^  V.^^)  /,  2. 

*  <^L1  "  '^L2>  /  ^• 
(R*)g2  =  3600. 


/  (12.  p.j'l 


(Re) 


,  =  3600.  p  -  V'  D  ,  /  (12.,  u  ,) 
rg2  '^g2  g2  e2  '  ''g2' 


1/2 


Call  FORCE  (D^.  D^,  V^,  V'^.  (S.^)j^,  (Re)^j.  (Re)^^, 

<*^*>rgr  <'^*'rg2’  Pgr  Pgr  «'■  ^i’  V 

F^  =  3600. 

<*c>1.2  ^  ’^'l2>  /  <^-  *’> 

If  (V  -  V  J  C157.  C160,  ^160 
g  L 

F.  =  -F. 

1  1 


F  = 
m 


(C142) 

(C143) 

(C144) 

(C145) 

(0146) 

(C147) 

(C148) 

(C149) 

(C150) 

(C151) 

(C152) 

{C153) 

(C154) 


m 


(C155) 

(C156) 

(C157) 

(C158) 


'/I 
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F»  =  -  F' 
m  m 

F  =  F  +  F. 
m  1 

F/(S.^)j3  - 

=  3600.  (w^/S.)^  V'^^/g. 


=  3600.  pi^  V.1^3  -L2. 

3600.  (w,  +  w.,)  V,  ,  +  g'  (F  +  F>  ) 


D.  ,/(12. 
cl'  'Ll 


m 


LZ 


3600.  (w  +  w  _) 

'±j  cz 


3600.  w  ,  V  -  g'  (F  +  F')  t 
V  - _ gl  gl  w 

g2  “  3600.  w 


g2 


rZ 


V  -  V 
g2  L2 


1/2 


V 


L2 


^3 

Pl 


(Re) 


0,2 


LZ 


(Re) 


rg2 


V. 

(0.  046)  3600. 


g' 


L  -J 


■v 


g2 


V  -  V 
r2  L2 


If  (Sense  Light  1)  C141,  C169 

CallGASP3  (T._(k),  O  ,  a  ,  <y,  t*.  (C/k)  ,  (C/k)  ,  W  , 

i2  nc  V  V  nc  nc 

^v'  P'  Pvi2'  Pgi2> 

If  ((Re)  ,  -  9000.)  C170,  C170,  C171 

rg2 

»G2  =  <''*>rg2‘^'  <^^2'^'  <Pg2/Pgi2>‘’' 

go  to  C172 

hG2  =  “-P”  <^*>rg2'’'®  <'’^>g2'^'  <Pg2/Pgi2''’' 


<^G  •  P'nc)2  =  '^02  <Pg2°'2\2>  '(S"'g2/<^-^>g2) 
If  ((Re)^-  -  12,500.)  C173,  C173,  C174 

LjC 


1/3 


(C159) 

(C160) 

(C161) 

(C162) 

(C163) 

(Cl  64) 

(C165) 

(C166) 

(C167) 

(C168) 


(C169) 


(C170) 

(C171) 

(C172) 


{C188) 
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If  (T  --T._)  C190.  C190,  C118 

g2  i2 

T.,  =T  -0.1 

i2  g2 

go  to  C118 


"2  = 


(w^/S.)  = 


\  (W2  ) 


'®i>12  =  <*c>12/'*c/®i> 


^12  '  =' 


/  u.  (S.),2  1 

2 

^°L2"°L1^ 

1 

M 

L\-2^V^  i 

i 

— 1  r/2 


^2  '  +  ^12 

“2=<*L*’"c2>/*g2 

temperature  ratio  =  T  -/T,  _ 

g2  L2 

phase  velocity  ratio  =  V  ^/Y.  _ 

gt.  Li2 

^2=®U  +  (Si'l2 

Q  =  ci  +  Q 

2  1  ^12 

percent  of  vapor  condensed  =  100.  w  -/w 

c2  VO 

R  _  =  1545.45/((1..X  ,)  W  +X  ,W  ) 
g2  v2  nc  v2  v 

Call  GASP  (T  -(K).  w  w  W  W  ,  J,  C  K  _) 
g2'  '  nc  v2  nc  v  pg2  g2 


^g2  =  <«o  ^2  ^2  Tg2(^)) 


1/2 


(C189) 

(CI90) 

(CI9I) 

(C192) 

(C193) 

(.C194) 

(C195) 


(G196) 

(€197) 

(C198) 

(C199) 

(C200) 

(C201) 

(C202) 

(C203) 

(C204) 

(C205) 


(C206) 


We  have  now  finished  calculations  "for  state  2.  If  we  satisfy  the 
following  tests,  the  calculation  is  compete;  if  not,  we  reset  several 
variables  to  establish  state  1  as  the  beginning  of  the  next  mixing  section 
interval  and  recalculate  state  2. 

If  (T  ,)  C234,  C234,  C207 

gZ  LiZ 


If  (w  rw  ,)  G234,  C234,  C208 

VO  c2 

(C207) 

If  (T  ,-T,  C234.  C234,  C209 

gZ  liZ  6 

(C208) 

^fgc  1  ^fgc2 

(C209) 

Svl  ^  %v2 

{C210) 

^pncl  ”  ^pnc2 

(C211) 

C  ,  =  C  , 

Pgl  Pg2 

(G212) 

^1-^2 

(C213) 

(C214) 

'"'vl  "  *v2 

(C215) 

T  =  T 

LI  L2 

(C216) 

T  =  T 
gl  g? 

(C217) 

Tgi(R)  =  T^2(R) 

(C218) 

^il  =  ^i2 

(C219) 

V  =  V 
gl  g2 

(C220) 

V  -Y 
.  LI  L2 

(C221) 

V  =  V 

LI  L2 

(C222) 

VI  =  V' 

gl  g2 

(C223) 

=  ^c2 

(C224) 

gl  g2 

(C225) 
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°L1  “ 

(C226) 

°gl  '  V 

(C227) 

(C228) 

®il  ■  ®i2 

(C229) 

”  °2 

(C230) 

(Re)gl  =  (Re)^, 

(C231) 

<*'*>rgl  =  <‘''>rg2 

(C232) 

(C233) 

go  to  C85 

read  in  data  for  new  case  (C234) 


Computational  Procedure  -  Analytical  Model  lib; 

The  computational  procedure  for  Analytical  Model  Ilb  is  essentially 
the  same  as  that  for  Analytical  Model  lla.  The  differences  which  are  listed 
below  are  the  result  of  simplifying  the  relative  velocity  between  phases  such 
that 

V  =  V  -V  s  V  +  V'  (G235) 

r  g  L.  g  Jj 

Thus,  all  equations  in  the  computational  procedure  for  Analytical  Model  Ila 

which  contain  the  relative  velocities,  V  or  V  ,  are  revised  to  use  the 

g  i-i 

single  relative  velocity,  V^. 

Delete  equations  C41  and  C42. 

Equation  G55  becomes 

(Re)^j  =  3600.  Ml,)  (C55) 

Equation  C56  becomes 

(R*>ngl  =  Pgl  ""el  /  "gP 


(C58) 
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Equation  C63  becomes 

f  s  3600.  (w  /S.),  V 
m  c  1  1 


rl 


/g' 


Delete  equations  C64  and  C65. 


Delete  equations  C115,  C116  and  C151. 
Equation  Cl 54  becomes 

PgZ  /  («. 

Call  FORCE  with  in  lieu  of  Y«  .  ♦ 

L  g 

Equation  C155  becomes 

F  =  3600.  (w  (V  -VJ/g' 
m  '  c  12  g  L  ® 

Delete  equations  Cl 56,  C159f  and  Cl 62. 

Equation  C163  becomes 
(Re)^2  '  3^“®-  Pl  ''tZ  °VZ  ' 

Equation  Cl 64  becomes 

3600.  (Wj^  +  w^j)  +  (g'F) 

^L2  ”  3600.  (w_  +  w  -) 

Li  cl 

Delete  equations  C167  and  C168. 

Delete  equations  C222  and  C223. 


(C63) 


(C154) 


(C155) 


(C163) 


(G164) 


Subroutine  Force  is  revised  slightly  to  account  for  this  change. 
{See  Appendix  G. ) 


57 


APPENDIX  D 

CONSTANT  AREA  MIXING  WITH  NONCONDENSABLE 
(ANALYTICAL  MODEL  ID) 


Computational  Procedure: 

We  begin  by  calculating  the  following  quantities  at  the  ejector  inlet 
state  (state  o)  and  at  the  entrance  to  the  iriixing  section  (state  1). 


I 

} 

i 


i 

f 


STATE  =  1. 


T_  =T  /(T  /T_  ) 

Lo  go  go  Lo 

Plo  =  Pgo  <PLo'Pgo> 

A  =w  /(w  +w  ) 
^  o  VO  VO  nc 


♦  1  =  ♦ 


w  ,  =  w 
vl  VO 


MV  SMV  +  W 
go  nc  VO 


w  ,  =  w 

gl  go 


T  (R)*  =  T  +  459.  69 
go  go 

T  (K)  =  5.  T  (R)/9. 
go  go 


X  .  = 


w  ./W 
vl  V 


vl  w  ,  w 

vl  nc 

W~  ^  W 

V  nc 


Call  GASP**  (T  (K),  w  .  w  , 
go  nc  VO 


W  .  W  ,  R',  C  .  K  ) 
nc  V  pgo  go 


(Dl) 

(P2) 

(D3) 

(P4) 

(D5) 

(D6) 

(D7) 

(D8) 

(D9) 

(DIO) 

(Dll) 

(D12) 

(D13) 


*  All  temperatures  are  in  degrees  F  except  as  indicated  by  (R)  or  (K)  which 
refer  to  the  Rankine  and  Kelvin  temperature  scales,  respectively. 

**  Fortran  Subroutines  are  identified  in  Appendix  G. 
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C  ,  =  C 
Pgl  Pgo 

K  =K 

gl  go 

R  ,  =  1545.45  /  ((1.  -  X  J  W  +  X  ,  W  ) 
gl  vl  nc  vl  V 


gl 


2.C,  JC  tT  (R) 

dff  iHE  O  ffO 


K  ,  R  ,  +  2.  C,  ^  JC  , 
1  gl  gl  dg  pgl 


-  459.69 


T  ,(R)  =  T  ,  +  459.69 
gl  gl 

Call  GASP  (T  ,(K).  w  .  w  W  ,  W  .  R*.  C  K  J 
gl  nc  vl  nc  V  pgl  gl 

‘‘^pgo  ^go''^)  -  ‘^pg,  ^gl!”»>] 

Si ' 


(D14) 

(D15) 

{D16) 


(D17) 


(D18) 

(D19) 


(D20) 

(D21) 

(D22) 


K  -1. 


Pl=Pgo^  "  ‘ 

*1^  =  3600.  PlVi^,Aj^j/144. 

CallSATT  (pj,  T^j. 

Call  CASPl  a  .  a  ,  a,  t*,  (c/k)  ,  {(/k)  ,  W  ,  W  , 

gl.  nc  V  '  'v  *  'nc  nc  v 

Pr  Si’  Sr  '^'r  '"vi-  '^ncr  '‘ei> 


(D23) 


(D24) 

(D25) 


Call  CASP2  {pj.  X^j,  W^,  R*.  T  (R).  T  (K). 


u  ,,  S'  ,  S'  ,  p  k  ,,  C  ,,  C  ,,  C  J 
^vl  nc  V  gl  ^gl  gl  pgl  pncl  pvl 


Si '  pgl  Si> 


(D26) 
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I>gl  =  (Agi  + 

(027) 

“  1  - 

(D28) 

CJ  =  &}  , 

o  1 

(D29) 

(D3b) 

Lj=0. 

(D31) 

S.,  =  0. 
ll 

(D32) 

Qj  =  0. 

(D33) 

phase  velocity  ratio  =  V  ,  /V,  , 

gl  L.1 

(D34) 

bulk  temperature  ratio  =  T  /T 

gl  Ll 

(D35) 

(Pr)  ,  =  C  ,  U  ,  /k  , 
gl  Pgl^  gl  gl 

(D36) 

Call  VISLIQ  (T^j. 

=  1  \l  -  \ 

(D37) 

**’''*L1  * 

{D38) 

°el  =  +  <=2  I’lIO 

(D39) 

(R*)gi  =  Pgi  '^gi  »^gl) 

(D40) 

'®=>gl=''g«''Pgl  °'l> 

(D4i) 

(Wc/S.)j  =  0. 

(D42) 

T.,  =  0.85  T  , 

ll  gl 

(D43) 

V  ,  =  0.75  V  , 
gl  rl 

(D44) 

Y\  =  0.25  V  , 

Ll  rl 

(D45) 

^il2  =  ®- 

(D46) 

Aqi2  =  0- 

(D47) 

E  =  0. 

(D48) 

60 


T.  j(R)  =  T.  j  +  459-  69 

Y. j(K)  =  5.  T. j(R)/9. 

CallSATP  (T.,,  p 

il  vil 

Call  SATHFG  (T. , .  p  . , .  h  , .  h.  ,) 
ll  Vll  gel  tel 

Call  yiSLIQ  (T.-  ,  u  ,) 

11  el 

^fgel  ^gcl  ^fel 

(Re)^j  =  3600.  ,  Mt  i) 


Ll  Ll' 


If((Rc)_  -  -  12,500.)  D53.  D53,  D54 

i-ll 

•'Ll  = 

go  to:T)55 

Hli  =  O  O”  (13.  k^/D^j)  »  (Pr)j^j®  '‘ 


(Re) 


rgl  =  Pgl  ^’gl  ?’.l/<‘^-  (‘gl) 


CaU  GASP3  (T  (K),  a  a  .  %  i*.  («/k)  .  (*/k)^^.  W  W  , 
1 1  ne  V  V  ne  ne  v 


If  URe) 


rgl 


P-,  p  u  .,) 

1  Vll  '^gll 

9000.)  D56,  D56,  D57 


■^Gl  =  (>'‘>rgl‘''  <^l/''gU>”'*^ 


go  to  D58 


0.8 


iK 


T  =  3600.  (w  /S.),  V«  ,/g' 
m  c  1  1  gl  ** 


V 


Ll 


> 


I  (Re) 


0.2 


Ll 


(Re) 


rgl 


gl  . .  2 


-I  1/2 


(0.046)  3600.  p 


L  «J 


(D49) 

(D50) 


(D51) 

(D52) 


(D53) 

(D54) 

.(D55) 


(D56) 


/  !  lO. 14 

<Mgi/Mgii) 

(D57) 

gl^<*’'''gl* 

1/3 

(D58) 

(D59) 

(D60) 
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v  =  V  -V 
gl  rl  LI 


ncl 


14. 696 


\  Pr^vi^i ' 


*^G1  ■  '  *’'nc*l^*’'ncl 


fgcl  pL'*il  *L1' 


~\l  ^^il‘^Ll^ 

Aqi  ^  .  (q^/S.)^ 

^ill  *  ^il2 
Aqi  1  -  Aq^2 

If  (E-50.)  D70.  D70,  D78 


2  2 
V  -V 
'  LI 


2.goJ 


-€^)  D78.  D78.  D71 


T  =  T 
il2  il 

Aqi2  =  Aqj 

E  —  E  +  1 . 


If  (E.-3.)  D74,  B75.  D75 

^il  "  '^142 
go  to  D76 


(D61) 

(D62) 

(D63) 

(D64) 

J  (D65) 

(D66) 

(D67) 

(D68) 

(D69) 

(D70) 

P71) 

(D72) 

(D73) 

(D74) 

(D75) 

(D76) 

(D77) 
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go  to  P49 


=  ^gl  -  ’■li 

Ui  =  (qL/Si)i/ATj 


(D78) 

(D79) 


At  this  point  we  have  calculated  all  variables  at  state  1.  We  now 
commence  our  calculations  for  the  first  etate  2. 


D  =  D 

g2  gl 

D  =D  . 
g  gl 

STATE  =  STATE  +  1. 


(D80) 

(D81) 

(D82> 


AT  =  c,  AT 
gig 

If  (AT  -€,)  D85.  D84,  D84 
g  2 

AT  =  € 
g  2 

T  ,  =  T  ,  -AT 
g2  gl  g 

If  (T  ,-T,  -€,)  D86,  D87,  D87 
g2  Lil  6 

^g2  -  "  ^6 

T  ,(R)  -  T  ,  +  459.69 

g2'  g2 

T  ,(K)  -  5.  T  ,/9. 
g2  g2 

Calls  ATP  (T^2«  Pv2^ 

P2  = 

If  (X  ,-X  ,)  D91,  D91,  D92 

V;1  v2 

T  =  T  -  1 

g2  g2 

go  to  D87 


(D83) 


(D84) 

(D85) 


(D86) 

(D87) 


(D88) 


(D89) 

(D90) 


(D91) 


V  =  V 
L2  LI 


(D92) 
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(D93) 

(D94) 

(D95) 

(D96) 

(D97) 

(D98) 

(D99) 

(DlOO) 

(DlOl) 

(D102) 

(D103) 

(D104) 

(D105) 

(D106) 

(D107) 

(D108) 

(D109) 


(DliO) 
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(w  ) ,  -  =  W  -  w 
c  12  gl  g2 


w  „  =  w  -  w  _ 
c2  VO  v2 


w 

If  (  -  €-)  D113.  D233,  D233 

w  3 

VO 

R  ,  =  1545.45/((1. -X  ,)  W  +  X  ,  W  ) 
g2  v2  nc  v2  v 

Call  GASP  (T  _(K),  w  ,  w  W  ,  W  ,  R',  C  K  ,) 
g2  nc  v2  nc  v  pg2  g2 


”  ,  -  (g  K  ,  R  ,  T  ,  (R)) 

'»o  g2  g2  g2  '  " 


1/2 


Call  GASPl  (T  _(K).  a  .  ff  ,  Of.  t*,  (c/k)  .  (c/k)  .  W  ,  W  . 

'  g2  nc  V  y  nc  nc  v 

Pr  *v2’  %2-  ®'2-  ''vZ'  "ncz-  V’ 

Call  OASP2  (p^,  R..  T  (R),  T  (K), 


S'  S'  ,  P  k  C  C  C  ,) 
nc  v  g2  ^g2  g2  pg2  pnc2  pv2 


'  Sg2  »‘gZ  ' 

'  Svl  ’^gl""  -  ^pv2 

(q  =  w  Ah 
V  12  v2  V 

Ah  =  C  ,  T  .  (R)  -  C  _  T  ,(R) 
nc  pncl  gl  pnc2  g2 

(q  )  =  w  Ah 

nc  12  nc  nc 

<S=)g2  = 

CallSATT  (p^, 

Call  SATP  (T.^.  p^.^) 

Call  SAtHFG  (T.^,  p^.2.  ^gc2’  ^fc2^ 

Call  SATLQV  (T._,  v  .) 

i2  c2 

Call  VISLIQ  (T.^  p.^^) 


Pc2  =  '•  ' 


c2 


(Dill) 


(D112) 


(P113) 

(D114) 


(D115) 
(DTI  6) 
(D117) 
(D118) 
{D119) 
(D120) 


(D121) 
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^fgc2  ^gc2  ^fc2 

AT  ^  =  (T  ,  +  T  _  -  T.  -  T.-)/2. 
gc  gl  g2  il  xZ” 

Ah  =  (C  .  +  C  J  AT  /2. 

VC  pvl  pv2  gc 

^fgc  ^  ^\cl 

Ap  =  144.  (p^-P2)/(PlJ) 


Ah  ^  =  C  . 
cL  pij 


T  +  T 
•’  i2 


/  il 
(  2- 


-  T 


L2 


) 


(q  )i,  =  W  (Ah  +  A,.  ) 

V  12  v2  V  Vg 

<V>1E  =  *nc 

'■Ic’u  =  <'*c>12  <\c  +  ^  ^^gL> 

<^L>12 

*^12  “  *‘’nc*12  '*'  *%*12  '''c'u  '*'  *^L*12 


’^L2  '  '^Ll  ■''C  ,  (w,  +  w  ,) 


pL 

AT  =  T  -  T 
2  g2  *L2 


cT 


(T  -T  -  T  -T 
at  -— g2  L2 
M  ,  T  -T, 


In 


(  gl  lM 
V  T  -T  i 

\  T  7  / 


g2  L2 

CaUVISLIQ(T^2* 


(D122) 

{D123) 

(D124) 

(D125) 

(D126) 

(D127) 

(D128) 

(D129) 

(D130) 

(P131) 

(D132) 

(D133) 

(D134) 

(D135) 

(D136) 

(D137) 

(0138) 


^^^^L2  "  %h  ^hZ^^L 


(D139) 


lifiSte-,--.  ^ 
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{w_/S,)  • 


c  »•  [  (w^/S.)j 


In 


\ 


-  1«. 


c2 


^c2^L2 


)/ 


3600. 


D.  .  =  (4.  A^/(»c,n 


1/2 


L2 


°e2  " 

V  ^  s  I  V  ,-V 


r2 


^2  = 


g2  ’L2 


‘®ic'l2  ^ 

2 

1  D,  ,-D,  , 
L.2  Lil 

2 

,  ‘'e2+*’i.1 

'*=2  2.  ' 

J 

1/2 


V  =  (V  -+V  _)/2. 

g  gl  g2 

(R*)g2  '  0g2\2°.2/<>"-  V’ 

(Re)^g2=  3600.,Pg^V.^^D^^/(12.Mgg) 

Call  FORCE  (Dj^,  D  Lj^.  V^.  V’^,  (S.j.)j2,  (R*)g2’ 

'rgr  <^*'rg2-  "gl-  V 


(Re)^^,.  {Re)_^,,  p^,,  p^,.  g'.  F,.  F  J 


F  =  3600.  (w  V  /g' 
m  cU2  g  ® 


(D140) 

(D141) 

(D142) 

(D143) 

(D144) 

(D145) 

(D146) 

(D147) 

(D148) 

{D149)' 

(D150) 

(D151) 

(D152) 

(D153) 

(D154) 


(D155) 
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F-^  =  3600.  + 


(D156) 


If  (V  -V^  )  D157,  D160,  D160 

fi  ^ 

T.  =  -F. 

1  1 

F  =  -F 
m  m 

F'  =  -F‘ 
m  m 

F  =  F  +  F. 
m  1 

^  =  I  ^/'®ic>12  I 

=  3600.  (w_^/S.)^  V^^/g. 

(Rc)i^^  ~  3600.  ^'l2  ^L.2^ 


(D157) 

(D158) 

(D159) 

(D160) 

(D161) 

{D162) 

(D163) 


P1+P2 


Vl2  = 


(PlALi-P2Ai^2)g'+3600.(Wj^+Wcj)Vi^j+g*(F+F'J.-^(Ai^j-A^2)g 


3600.  (w_  +w  ,) 
cz> 


(D164) 


V  S  V  -V 
rZ  gZ  LZ 


(Re)j 


V  ^  + 

g2 


(0. 046)  3600. 


V'  =  V  _  -  V'  , 
gZ  rZ  LZ 


P2  =  Pi  + 


3600- 


If  (Sense  Light  1)  D107,  D169 


(DI65) 


(DI66) 


(D167) 


(D168) 


CallGASP3  (T.2(K).  «.  t*.  (c/k)^.  (c/k)^^,  W^. 

^Z’  ^viZ’  ^giZ^ 

If  ((Re)  ,-9000.  )  JD170.  D170,  D171 

rgZ 


(D169) 
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go  to  D172 


“gE  =  ®  V'^’eE*  ‘**'>,2^^^  ‘  V'^giz’”' 

<‘'g  •  *"nc>2  -  “gZ  V'^’Z  V  «S'>gz/'*’V‘'" 

K  ((RcK  --12.500.)  D173.  D173.  D174 

XjZ 

=  1.86  (I2.kj^/Dj^)  (Re),  ,*'*  (Pe),  ,*'*  (m,  ,/M,/  ** 


'L2 


'L2 


‘L2''‘c2' 


go  to  D175 
h 


L2  '  ®  ®“  ““’lz"'®  <*’^>L2‘’  ‘‘ 


P* 


^v2^2"^vj2 


nc2 


14. 696 


In 


Pa-^vzPz 


‘‘gZ  '  ‘•^G  •  *”„c'z'*”ncz 
<’*c'®i»Z=»‘Gz'*vzPz-Pv«'/‘*<”‘ 


<''g'®i>z“‘'GZ<V*^i2>^<*c'®i'z{’ 


V.  -V 

^fgc2'‘’SL^^i2‘^L2^^  *2g  J 


”  ^L2  ^^i2'^L2^ 
Aq,  =  (q^/Sp^  -  (q^/S.)^ 

^i21  “  ^i22 

Aq^i  =  Aq^^ 

If  (E-50.)  D183,  D183,  D191 


Ii(l 

'^iz- 

T 

S  T 

122 

12 

^^22 

^Aq^ 

-fp  D191.  D184,  D184 


(D170) 

(D171) 

(D172) 

(D173) 

(D174) 

(D175) 

(D176) 

•(D177) 

\ 

0 

-j|(D178) 

(D179) 

(D180) 

(D181) 

(D182) 

pi  83) 
pi  84) 


pi  85) 
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(D186) 

(D187) 

(D188) 

(D189) 

(D190) 

(Din) 

(D192) 

(D193) 

(D194) 

(D195) 

(D196) 

(D197) 

(D198) 

(Di99) 

(D200) 

(D201) 

(D202) 

(D203) 
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We  have  now  finished  calculations  for  state  2.  If  we  satisfy  the 
following  tests,  the  calculation  is  complete;  if  hot,  we  reset  several 
variables  to  establish  state  1  as  the  beginning  of  the  next  mixing  section 
interval  and  recalculate  state  2. 


if(T  D233.  D233.  D204 

gZ  liZ 

If  (w  -w  J  D233.  D233.  D205 

VO  c2 

(D204) 

If(T  J  D233.  D233.  D206 

gZ  L2  6 

(D205) 

^fgcl  *  **fgcE 

(D206) 

Svl  *  SvZ 

(D207) 

^pncl  ”  ^pnc2 

(D208) 

Pgl  Pg2 

(D209) 

Pi  *  P2 

(D210) 

(w  IS.).  =  (w  IS.). 
c  rl  '  c  X  2 

(D211) 

^vl  "  %2 

(D212) 

^cl  ~^c2 

(D213) 

w  .  =w 
gl  g2 

{D214) 

T  s  T 

LI  L2 

(D215) 

*  V 

(D216) 

Tgi(R)  = 

(D217) 

^il  ^  ^i2 

(D218) 

V  =  V 
gl  g2 

(D219) 

V  =  V 

LI  L2 

(D220) 

V  ,  =  v« 
gl  g2 

(D221) 
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V  s=  V‘ 

LI  L2 

(D222) 

D  s  D 

LI  L2 

(D223\ 

L,  =  L, 

1  2 

(D224) 

Q  Q 

1  2 

(D225) 

(Re)gi  =  (Re)^ 

(D226) 

(D227) 

^L1  ”  ‘^L2 

(D228) 

A  =  A  , 
gl  g2 

(D22,9) 

S  ,  s  S 
il  i2 

(D230) 

^vl  *  ^y2 

(D231) 

(D232) 

go  to  D82 

read  in  data  for  next  case  (D233) 
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APPENDIX  E 

DERIVATION  OF  SOME  BASIC  CONSERVATION  EQUATIONS 


Continuity  Eqtiationa : 


In  our  mixing  section  analyses  we  assume  that  two  flows  enter  the 
ejector.  The  secondary  or  liquid  stream,  w  ,  and  the  primary  or  gas 

Lj 

stream,  w  ,  which  may  consist  of  a  condensable  vapor,  w  ,  or  a  mix- 
go  '  VO 

ture  of  condensable  vapor,  w  ,  and  k  noncondensable  gas,  w  .  The 

VO  *  nc 

noncondensable  flow  rate,  w  ,  remains  constant-throughout  the  mixing 

nc 

process,  whereas  the  vapor  flow  rate  generally  decreases  because  of 
condensation.  That  portion  of  the  vapor  which  condenses  is  defined  as 
w  .  The  total  amount  of  condensed  vapor  at  state  1  is  defined  as  w  .  and 

C  C  X 

the  vapor  which  condenses  between  any  two  states  1  and  2,  as  (w  )... 

c  Ic 

Thus,  at  state  1  we  define  ^e  liquid  flow 

LI  L  cl 

and  the  gas  flow 

w  ,  s  w  ,  +  w 
gl  vl  nc 

At  state  2  the  liquid  flow  becomes 

*L2'*Ll+<'"c>12 
and  the  gas  flow 

w  _  s  w  -  -  (w  ) .  „ 
g2  gl  c'12 

Also,  the  total  condensed  vapor  at  state  2  is 


(El) 


(E2) 


(E3) 


(E4) 


c2  cl  c  12 

For  the  present  mixing  section  analyses,  multiple  liquid  jets  may 
be  assumed.  In  this  case  we  distribute  the  liquid  flow  evenly  between  the 
jets  and  assume  that  condensation  occurs  equally  on  each  jet  at  any  given 
state,  If  we  let  c^  equal  the  number  of  liquid  jets,  then  the  liquid  jet 
diameter  becomes  at  state  1 


(E5) 
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4,  A 


1/2 


= 


1^1 


Jfc 

"^2 


(E6) 


where  the  total  flew  area  lor  liquid,  is 

A  _  144.  / 

U  =  36X,0 

^he  flow  area  for  the  gas  mixture  is 


(E7) 


A  ,  = 


144. 


w 


ii 


gl  3600.  1  p^jV^, 

.and  the  diameter  of  the  mixing  section  becomes 


(E8) 


D 


gl  I  tr 

Thus,  for  any  given  state  in  the  mixing  section,  we  define  the  total  flow 
area  as  being 


(E9) 


A  =  A,  T  A 
E  g 


(ElO): 


Momentum  Equations,  Constant  Pressure  Mixing; 

With  reference  to  Figure  27a,  the  conservation  of  linear  momentum 
for  the  liquid  jet  control  volume  yields 


?l'^Ll‘^ 


3600.w^^V^l 

g' 


+  F.+F  +F'  =p,A 


m 


m  2  L2 


3600. 

+ - p - + 


Pl+Pj 

2. 


(Eli) 


where  we  assume  a  linear  pressure  variation  in  the  mixing  section  and  have 
neglected  any  angle  effects  on  the  force  along  the  sides  of  the  control  volume. 

From  Figure  2  7b  the  conservation  of  linear  momentum  for  the  gas 
mixture  control  volume  yields 


3600.W  V  p^+P2 

PlAgi+  g.  +  2.  ^^Lr^L2^''P2\2‘*’ 


3600.W  V  , 

- - +  F  +F.+F  + 

g'  w  1  m 


(P1+P2)  (Aj-A2)/2. 


(E12) 
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Referring  to  Figure  28,  the  total  relative  velocity  between  the  bulk 
gas  mixture  and  the  bulk  liquid  stream,  v^,  i& 

V  sV  -V, 
r  g  L 

also. 


(E13) 


V  -V*  (E14) 

r  g  Li 

where  is  the  relative  velocity  between  the  bulk  gas  mixture  and  the 

liquid  jet  surface,  or  phase  interface,  and  V'  is  the  relative  velocity 

between  the  phaf.&  interface  and  the  bulk  liquid  jet.  To  calculate  the  rela> 

tive  velocities  V  and  V  at  any  state,  we  proceed  as  follows.  From 

g  L> 

shear  stress  at  the  phase  interface  we  get 

T  +r  (E15) 

g  m  L 

where  T  is  shear  stress  due  to  momentum  exchange  across  the  interface 
xn 

and  T  and  T,  are  shear  stresses  due  to  friction  at  the  interface, 
g  L 

Since 


3600 


■Ct) 


g 


T  s 
m 


g' 


(El  6) 


and  from  Reference  20 


(f.)  py'  3600.^ 

T  _ 

g  2.g- 


(E17) 


2.g' 


(E18) 


where  the  Fanning  friction  factors,  (f.)  and  are  calculated  from 
equation  E33  (streamline  flow)  and  equation  E 35  (turbulent  flow),  and 
the  Reynolds  nujmbers  are 
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.  3600.  / 


p  V‘  D 

^  g 


(E19) 


to  calculate  (f.)  and 
1  g 


,  3600.  f  \ 

{-H— j 


(E20) 


to  calculate  (f.)  ,  we  may  solve  equations  E15-E20  for  V  with  the  result 

1  Li  Li 


5!lP) 

(f;),  \pJ 


i'L  ''"L 

and  from  equation  E14 


2 

V  + 
g 


2-gT 


1/2 


m 


(f.>^p^  3600. 


(E21) 


v*  «  V  -V\  (E22) 

g  r  L 

Having  determined  the  relative  velocities  V'  and  V  ,  we  may  now  calculate 

^  g 

the  following  resultant  shear  forces  from, 

momentum  exchange  from  the  bulk  gas  mixture  to  ti^te  phase  interface 

F  =  3600.  (w  V'  /g'  (E23) 

m  c  12  g  ® 

mpmentum  exchange  from  the  phase  interface  to  the  liquid  jet 

F'  =  3600.  (w  V'  /g'  {E24) 

m  c  12  Li 

wall  friction 


f  p  V  ^  3600.  ^ 
yr  g  g _ 

w  >,\  2.g' 

and  phase  interface  friction 


^^wh2 


F.  = 
1 


'  f.p  V  ^  3600.^ 
2.g’ 


<"i>12 


\  / 

Equations  E25  and  E26  are  approximations  of  the  friction  forces  in  the 
x-direction  given  by 


(E25) 


(E26) 


'  Tt- . * 
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and 


M  T  D  dL  C08  P 
w  g 


ff  T.  D,  dL  cos  Y 
1  L 


(E27) 


(£28) 


where  we  have  assumed  steady,  one -dimensional'  flow,  and  uniform  pres¬ 
sure  across  sections  1  and  2,  For  purposes  of  the  present  calculations 
we  also  assume  that  the  half-angles,  fi  and  Y,  are  small.  Thus,  we  may 
rewrite  equations  £27  and  £28  as 


^  (S 

w  w  w  12 


(E29) 


and 


F.  (S.),  _ 

1  1  1  12 

In  accordance  with  McAdams  (Reference  20)  we  define  the  shear  stress  as 


(£30) 


f  p  3600,^ 

2.g' 

and  we  estimate  friction  factors  (Reference  20)  as  followsf : 

streamline  flow:  (Re  1000) 

f  =  16.  /(Re) 
w  g 

f  =  16. /(Re) 

1  rg 


(£31) 


(£32) 

(£33) 


turbulent  flow:  (Re  1000.) 

f^s  0.0014+  /o.  125/(Re)  (E34) 

f.  =0.0014+  ^0.  125/(Re)^^°*^^j  (E35) 

The  mixing  section  wall  friction  factor  is  based  on  the  bulk  gas  mixture 
Reynolds  number 


p  V  D  3600. 

(Re)  =^V . - 

g  12.  Mg 


(£36) 
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In  the  constant  pressure  mixing  section  ah^yses  for  Analytical 
Models  I  and  Ila,  where  from  constant  pressure 

P"Pl=P2 

we  can  solve  equation  Ell  for  the  liquid  jet  velocity 


3600.  +  g.  (F^  +  F.  +  F'_^) 


L2 


3600.  w 


L2 


and  combining  equations  ElO  and  E12  we  can  solve  for  the  bulk  gas  mix 
ture  velocity 


3600.  w  V  -  g'  (F_  +  F'  +  F  ) 
_ g  1  =g  1  m  1  w 


3600.  w 


g2 


In  the  constant  pressure  mixing  section  analysis  for  Analytical 
Model  Ilb  where  we  have  simplified  the  relative  velocity  to  the  expression 
in  equation  E13,  equation  E3.8  becomes 


3600.  *  g.  (F  +  F.) 


L2 


360,0.  w 


L2 


and  the  shear  forces,  F  and  F.,  are  evaluated  as  follows: 

m  i 

F  =  3600.  (w  k  V  /g' 
m  c  r 


F.  = 
1 


f.  p  V  3600. 
1  r _ 

2.  g' 


where  the  average  friction  factor,  f^,  is  evaluated  on  the  basis  of  the 
Reynolds  number,  (Re)^^,  given  by  equation  E36  and  redefined  as 


(Re) 


p  V  D  3600. 
_  g  r  e 


rg 


12.  p 


g 


(E37) 

(E38) 

(E39) 

(E40) 

(E41) 

(E42) 


(E43) 
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Momentam  Equation* ,  Coaatant  Area  Mixing: 

In  the  mixing  section  analysis  for  Analytical  Model  Ill  where  we 
assurne  that  the  mixing  section  remains  at  constant  area 

we  solve  equation  £11  for  the  liquid  jet  velocity 

-I  / 

and  the  bulk  gas  mixture  velocity  from  equation  £  1 1 


V 


^  144. 
g  "  3600. 


) 


and  tba  total  static  pressure  in  the  mixing  section  from  equations  £12  and 
£44 


P2  =  Pi  "■ 


3600.  {w„,V^,-w  V  )-g>  (F  +r+F  ). 
gf  ■  gl  8^  8^  m  1  w 


g' 


( 


'^Ll'^^LZ 


+  A 


g2 


) 


The  resultant  shear  forces,  F  ,  F'  ,  F  and  F,  arc  evaluated  from 

m  m  w  1 

equations  E23,  D24,  E25  and  E26  respectively. 


Energy  Equations; 

The  steady  flow  energy  equation  for  the  control  volume  represented 
by  Figure  29  may  be  written  as 


w 


gi 


/h  +w,  ,  / 

[  gl  Z-g^J/  LI  [ 


h 

LI  2.g  J 


)  g2  i  g2  2.g^j; 


V, 


w 


L2 


L2 


L2  2.g  J 
*o 


(£44) 


(E45) 


(£46) 


(E47) 


(E48) 
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Since  from  equations  £4  and  £3. 

(E49) 

(E3) 

by  combining  equations  E48,  E49.  and  E3  we  rewrite  the  energy  equation 
as 


and 


w  1  li  — h  4-  —  fit 

gZ  ^  gl  gZ*  z 


2  2 
V  .  -V 


17^  )•'-.■»( 


^^Vl2  rvl'^L2‘^^ 


2  2 
’^L2  \ 

2g  J  / 


2  2 
V.  .  -V. 


w 


Ll 


V  ^  Ll  '  L2  \  ^ 

(*'Lr^L2^  2.g  J  j  ‘ 


Assuming  that 


(E50) 


l"^T  ^  C  T  (-T  l"T_  _)  + 


144. 


(Pi'Po) 


Ll  L2  pL  '*L1  L2'  p  J  '*^1  ^2 


(E51) 


and  substituting  equation  E51  into  E50  and  solving  for  T  ,  our  energy 

LfZ 

equation  now  becomes 


T.  =  T.  .  + 


1. 


L2  Ll  w,  ,C  ^ 
Ll  pL 


w  /h  h  ) 

g2  \  gl  g2  2.g  J  / 


+  <*c>12 


+ 


2.g 


J 

o 


<PrP2^U 

{E52) 


Now  for  the  mixing  section  interval  in  Figure  29  we  define  the  following 


individual  heat  transfer  rates. 
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(q  the  heat  transfer  rate  from  die  ncmcondensable 

nc  12 


(2  2 

V  -V  - 

C  ,T  ,(R)-C  ,T  .(R)+-8i - 

pncl  pnc2  g2  2.g  J 


(q  the  heat  transfer  rate  from  the  vapor  which  does  not  condense 

V  12 


(E53) 


(E54) 


(q  )  ,  the  heat  transfer  rate  from  the  vapor  which  condenses 

C  J.  M 


<Vl2  =  <*c>12 


(C  ,+C  ,){T  ,+T  ,-T.  -T.,)+h.  (T.,)  +  h-  (T.,) 

^  pvl  pv2  gl  g2  il  i2  fg^  il  fg  i2 


2  2  \ 

,  T.-+T._  .  V  /-V.-  \ 

[  z.  ^vz)  *  z.g^J  I 


(q  )  ,  the  heat  transfer  rate  from  the  liquid  jet 

±j  IZ 

(2  2 
V  -V 

^L1  LZ  ^  144.  ,  , 

2.goJ  p^j  ^^r^2^ 


and  the  total  heat  transfer  rate  for  the  interval  1-2  as 


(ESS) 


(E56) 


^12  ‘  ^%chz  ^  ^Vl2  ^Vl2  ^^L^12 


(E57) 


Substituting  equation  E57  into  equation  E52  yields 


T  —  'p  j.  *** 

L2  "  LI  w,  ,C  , 
LI  pL 


(E58) 


For  the  constant  pressure  mixing  section  analyses,  Analytical  Models 
I,  Ila,  and  lib  the  pressure  term  in  equation  E56  drops  out  to  yield 
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‘’i’lZ 


For  Analytical  Model  I  where  w^^  >  0. ,  equation  £53  is  deleted  and 
equation  £57  becomes 


^12  "  ^Vl2  ^  ^  ^*^L^12 


(£59) 


(£60) 


\ 

I 

I 
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APPENDIX  F 

SOURCES  OF  DATA  AND  CORRELATIONS 

Sources  of  data  for  property  calculations  and  of  correlations  for  heat 
and  mass  transfer  calculations  used  in  the  mixing  section  analyses  are  given 
below^ 

Property  Calculation  Parameters; 

A.  Water  Vapor 

1.  Lehnard- Jones  force  constants 

a.  molecular  potent^l  well  depth,  (C/k)  ,  356.  K 

V  • 

b.  molecular  collision  diameter,  0,2.  649  A 

V 

Table  8-2,  p.  270,  Reference  24. 

c.  measure  of  molecular  polarity,  t*,  1.2 
Table  3.  10-1,  p.  214,  Reference  23. 

2.  Sutherland  constant.  S'  ,  559.  74  K 

V 

equation  7-19,  Reference  24. 

B.  Carbon  Dioxide  Gas 

1.  Lennard- Jones  force  constants 

a.  molecular  potential  well  depth,  {</k)  ,  213.  K 

nc 

b.  molecular  collision  diameter,  O  ,  3.  996  A 

nc 

Table  I- A,  p.  Mil,  Reference  23. 

-23 

c.  molecular  polarizability.  Of,  0.265  x  10 
Table  13.2  -  3,  p.  950,  Reference  23. 

2.  Sutherland  constant,  S'  ,  233.  K 

nc 

equation  7-19,  Reference  24  (revised  to  agree  with  available 
vis  CO'S  ity  data)  . 
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Fanning  Friction  Factors; 

Mixing  section  wall  and  phase  interface  friction  are  both  estimated 
from  correlations  for  fluid  friction  in  pipes  based  on  isothermal  fluid 
flow. 

A.  Streamline  flow,  Re  ^  1000. 

f  =  16.  /Re  (FI) 

equation  3. 46,  Reference  19 
equation  6 -5a,  Reference  20 

B.  Turbulent  flow.  Re  ^1000.  (Smooth  wall) 

i  ‘  0.  0014  +  (F2) 

Re 

equation  3. 47e,  Reference  19 
equation  6.  8,  Reference  20 

For  purposes  of  the  mixing  section  calculations,  a  smooth  transition 
front  streamline  to  turbulent  flow  was  assximed  at  Re  =  1000; 


Heat  Transfer  Correlations; 

A.  Annular  gas  stream 

1.  Streamline  flow,  (Re)  S  9000. 

rg 


,  12. k  V 

(V) 


h^  =  1.86 

Li 


s  1  (Re)  (Pr) 


rg 


g 


equation  6.  1,  Reference  19 
equation  9-28a,  Reference  20 


2.  Turbulent  flow  (Re)  >  9000. 

rg 


=  0.027 


(Re) 


0.8 


1/3 


rg 


(Pr) 


0.  14 


(F3) 


(F4) 
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equation  6.  3,  Reference  19 


B.  Liquid  Jet 

1.  Streamline  flow,  (Re)  ^12, 500. 

Xj 


equation  6. 1,  Reference  19 
equation  9-28a,  Reference -20 


(F6) 


2.  Turbulent  flow,  (R^)^^  12, '500. 

Ij 

h^.  0.023  {^)  (Re)J-«(Pr)“-'‘ 
equation  9- 10a,  Reference  20 


(F7) 


It  has  been  assumed  for  these  calculations  that  the  heat  transfer 
coefficients,  h_  and  h^  ,  for  streamline  flow  are  independent  of  the  term 
L/D  appearing  in  the  correlations  given.>by  the  references.  Also,  to  suit 
the  calculations,  the  transition  points  from  streamline  to  turbulent  flow 
are  taken  as  (Re)^^  =  9000.  for  the  gas  stream  and  (Re)^  =  12,500.  for 
the  liquid  jet. 


Mass  Transfer  Correlation; 


(F8) 


equation  13.29,  Reference  19. 
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APPENDIX  G 

IDENTIFICATION  OF  SUBROUTINES  USED  IN 
COMPUTATIONAL  PROCEDURES 


Several  Fortran  subroutines  used  in  the  computer  calculations  and  which 
appear  in  the  computational  procedures  of  Appendices  B ,  C  and  D  are  des  - 
cribed  briefly  below. 

SATP; 

Description  -  steam  and  water  property  calculation  to  determine 
saturation  pressure  corresponding  to  a  given  temperature. 

Input  -  temperature,  T  F 

Output  -  pressure,  p  psia 

Limits  >  pressure  is  given  by  equations  11  and  12  of  the  reference. 
Reference  -  Reference  21. 


SATT: 

Description  -  steain  and  water  property  calculation  to  determine 

saturation  temperature  corresponding  to  a  given  pressure. 

Input  -  pressure,  p  psia 

estimated  temperature,  T  F 

Output  -  temperature,  T  F 

6  at 

Limits  -  trial  and  error  calculation  with  pressure  given  by 
equations  11  and  12  of  the  reference.  Specified  accuracy  for  solution 


derived  from  two  consecutive  calculations  is 


P2.P1 


<0.0001. 


Reference  -  Reference  21. 


VAPTP: 

Description  -  steam  and  water  property  calculation  to  determine 
enthalpy,  entropy,  and  specific  volume  of  vapor  at  a  given  temperature 
and  pressure. 
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Input  -  temperature,  T  F 

pressure,  p  psia 

Output  -  enthalpy,  Btu/lbm 

entropy,  s  Btu/lbin-F 

3 

specific  volume,  v  ft  /Ibm 


g 

Limits  specific  volume,  enthalpy,  and  entropy  are  given  by 
equations  13,  14  and  16A.  of  the  reference. 

Reference  -  Reference  21. 


SATLQ: 


Description  >  steam  and  water  property  calculation  to  determine 
enthalpy,  entropy,  and  specific  voliune  of  saturated  liquid  at  a  given 
temperature. 

Input  -  temperature,  T  F 

Output  -  enthalpy,  h^  Btu/lbm 

entropy,  Btu/lbm-F 


specific  volume,  v 


ft  /Ibm 


Limits  -  enthalpy  and  entropy  are  polynomial  approximations  based 
on  the  tabulated  values  of  the  reference  between  the  temperature  limits 
50-350  F.  Specific  volume  is  given  by  equation  18  of  the  reference. 
Reference  -  Reference  21. 


SATHFG; 


Description  -  steam  and  water  property  calculation  to  determine 
saturated  vapor  enthalpy  and  saturated  liquid  enthalpy  at  a  given  saturated 
temperature  and  pressure  condition. 

Input  -  temperature,  T  F 

pressure,  p  psia 


Output  -  vapor  enthalpy,  h 


g 


liquid  enthalpy,  h^ 


Btu/lbm 

Btu/lbm 
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Limits  -  liquid  enthalpy  is  a  polynomial  approximation  based  on  the 
tabular  values  of  the  reference  between  the  temperature  limits  50-350  F. 
Vapor  enthalpy  is  given  by  equation  14  of  the  reference. 

Reference  -  Reference  21. 

SATLQV: 

Description  -  steam  and  water  property  calculation  to  determine 
specific  volume  of  saturated  liquid  at  a  given  temperature. 

Input  -  temperature,  T  F 

3 

CXitput  -  specific  volume,  v^  ft  /Ibm 

Limits  -  specific  volume  is  given  by  equation  18  of  the  reference. 

Reference  -  Reference  21. 

SUPHTl; 

Description  -  steam  and  water  property  calculation  to  determine  a 
superheated  vapor  temperature  at  a  given  pressure  and  vapor  entropy. 
Inpv.t  -  temperature,  IJ?  (estimated)  F 

pressure,  p  psia 

entropy,  s  Btu/lbm-F 

Output  -  temperature,  T  F 

Limits  -  trial  and  error  calculation  with  entropy  given  by  equation 

16A  of  reference.  Specified  accuracy  for  solution  derived  from  two 

consecutive  calculation  is  “®  i 

— — -  <  0.0001. 

Reference  -  Reference  21. 

SUPHT2; 

Description  -  steam  and  water  property  calculation  to  determine  a. 
superheated  vapor  temperature  at  a  given  pressure  and  vapor  enthalpy. 
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Input  -  temperature,  T  (estimated) 
pressure,  p 


enthalpy,  h 


g 


psia 
Btu/lbm 

Output  -  temperature,  T  F 

Limits  -  trial  and  error  calculation  with  enthalpy  given  by  equation  14 
of  reference.  Specified  accuracy  for  solution  derived  from  two  consecutive 


calculations  is 


h^-hi 

"i 


<0.0001. 


Reference  -  Reference  21. 


VISLIQ; 

Description  -  calculation  to  determine  viscosity  of  water  at  a  given 
temperature. 

Input  -  temperature,  T  F 

Output,  viscosity,  U-  Ibm/hr-ft 

La 

Limits  -  tabular  values  from  32  to.  212  F  and  equation  are  given  on 
p.  374  of  reference. 

Reference  -  Reference  22. 


VAP: 


Description  -  calculation  to  determine  density  viscosity,  ratio  of  specific 
heats,  and  specific  heat  at  constant  pressure  of  water  vapor  at  a  given  temp- 
eratiire  and  pressure. 

Input  -  temperature,  T  F 

temperature,  T  K 

pressure,  p  psia 


molecular  potential  well  depth, 

(f/k) 

V 

molecular  weight, 

molecular  collision  dia¬ 
meter,  O’ 

V 

universal  gas  constant,  R' 


K 

Ibm/lbm-mole 


Btu/lbm-mole-R 
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Output  -  density,  Ibm/ft 

viscosity,  Ibm/hr-ft 

ratio  of  specific  heats ,  K 

g 

specific  heat  at  constant 

pressure,  S  Btu/lbm-R 

Pg 

Limits  -  Specific  volume  is  given  by  equation  13  of  Reference  21. 
Specific  heat  at  constant  pressure  is  given  by  Table  170,  p.  220  of  Ref¬ 
erence  22  for  the  temperature  range  300-2500  K.  Viscosity  is  given  by 
equation  8.2-18  of  Reference  23.  Collision  integrals  are  approximated 
from  Table  I-M  of  Reference  23. 

References  -  specific  volume,  Reference  21 

specific  heat  at  constant  pressure,  Reference  22 
viscosity.  Reference  23. 

GASP; 

Description  -  calculation  to  determine  specific  heat  at  constant  pres¬ 
sure  and  ratio  of  specific  heats  of  a  carbon  dioxide  gas  and  water  vapor 
mixture  at  a  given  temperature  and  known  flow  rates. 

Input  -  temperature,  T  K 

noncondensable  flow  rate,  w  Ibm/hr 

nc 

condensable  flow  rate,  w  Ibm/hr 

V 

noncondensable  molecular 

weight,  W  Ibm/lbm-mole 

*  nc 

condensable  molecular 

weight,  Ibm/lbm-mole 

universal  gas  constant,  R'  But/lbm-mole-R 

Output  -  specific  heat  at  constant  pressure  of  the  gas 

mixture,  C  Btu/lbm-R 

pg 

ratio  of  specific  heats  of  the 
gas  mixture, 

Limits  -  specific  heat  at  constant  pressure  of  carbon  dioxide  gas  and 
water  vapor  given  by  Table  170,  p.  220  of  the  reference.  Temperature 
range  for  gas  is  273-1200  K.  and  for  the  vapor  is  300-2500  K. 

Reference  -  Reference  22. 
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GASPl: 

Description  -  gas  and  vapor  mixture  property  calculation  to  determine 
gas  mixture  molecular  collision  diameter,  binary  diffusion  coefficient, 
condensable  viscosity,  noncondensable  yjscosityi  and  gas  mixture  viscosity 
for  a  given  temperature  and  condensable  mole  fraction. 

Input  -  temperature,  T  .  K 

none ondens able  molecular.. c pllis ion 
diameter,  O 

nc 

condensable  molecular  collision 
diameter^  O 

V 

property  calculation  parameter,  a 

property  calculation  parameter,  t* 

condensable  molecular  potential 
well  depth, 

noncondensable  molecular  potential 
well  depth,  (€/k) 

nc 

Aoncondensable  molecular  wfeight,  W 

®  nc 

condensable  molecular  weight, 

pressure,  p 

condensable  mole  fraction,  X 

V 

Output  -  gas  mixture  molecular  collision 
diameter,  <T 

g 

binary  diffusion  coefficient,  D' 

condensable  viscosity, 

noncondensable  viscosity,  u 

^  ^nc 

gas  mixture  viscosity,  n 

g 

Limits  -  individual  gas  mixture  component  viscosities  are  given  by 
equation  8.2-18  of  Reference  23.  Gas  mixture  viscosity  is  given  by 
equation  8.2-30  of  Reference  23.  The  binary  diffusion  coefficient  is 
given  by  equation  8.2-44  of  Reference  23.  Molecular  collision  dia¬ 
meter  and  molecular  potential  well  depth  of  the  mixture  are  given  by 
equations  8.6-3,  8.6-4  and  8.6-5  of  Reference  23.  Collision  integrals 
are  approximated  from  Table  I-M  of  Reference  23. 

Reference  -  Reference  23. 


A 

A 

K 

K 

Ibm/lbm-mole 

Ibm/lbm-mole 

psia 

A 

ft^  /hr 
Ibm/hr-ft 
Ibm/hr-ft 
Ibm/hr-ft 
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GASP2: 


Description  -  carbon  dioxide  gas  and  water  vapor  mixture  property 
calculation  to  determine  density,  thermal  conductivity  and  specific  heat 
at  constant  pressure  for  the  gas  mixture,  noncondensable  specific  heat 
at  constant  pressure,  and  condensable  specific  heat  at  constant  pressure 
for  a  given  temperature,  pressure,  and  condensable  mole  fraction. 

Input  -  pressure,  p  psia 


condensable  mole  fraction,  X 


noncondensable  molecxilar  weight,  W  Ibm/lbm-mole 

®  nc 


condensable  molecular  weight, 
universal  gas  constant,  R‘ 
temperature,  T 
temperature,  T 
noncondensable  viscosity, 
condensable  viscosity,  n 


Ibm/lbm-mole 

Btu/lbm-moie-R 


Ibm/hr-ft 


Ibm/hr-ft 


noncondensable  Sutherland  constant, S'  K 

nc 

condensable  Sutherland  constant.  S'  K 


gas  mixture  molecular  collision 
diameter,  O 

g 

Output  -  gas  mixture  density,  p 


Ibm/ft' 


gas  mixture  thermal  conductivity, k  Btu/hr-ft-F 

gas  mixture  specific  heat  at  constant 

pressure,  C  Btu/lbm-R 

Pg 

noncondensable  specific  heat  at 

constant  pressure,  C  Btu/lbm-R 

pnc 

condensable  specific  heat  at  constant 

pressure,  C  Btu/lbm-R 

pv 

Limits  -  specific  heat  at  constant  pressure  of  carbon  dioxide  gaS'  and 


Btu/lbm-R 


Btu/lbm-R 


water  vapor  given  by  Table  170,  p.  220  of  Reference  23.  Temperature' 
range?  for  gas  is  273-1200  K,  and  for  the  vapor  is  300-2500  K.  Thermal 
conductivity  of  binary  mixture  is  given  by  equation  7-17  of  Reference  24. 
References  -  specific  heat  at  constant  pressure.  Reference  23. 

'as  mixture  thermal  conductivity,  Reference  24. 
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GASP3: 

Description  -  gas  and  vapor  mixture  property  cal ctilation  to  determine 
viscosity  of  mixture  at  a  given  temperature  and  condensable  mole  fraction. 

Input  -  temperature,  T  K 

n^ncondensable  molecular  collision 

diameter,  ^  A 

~  ~  nc 

condensable  molecular  collision 

diameter,  Q  A 

y  " 

property  calculation  pafanieter,  a 

property  calculation  parameter,  -t* 

condensable  molecvilar  potential  well 
depth,  (C/k)^  K 

noncondensable  molecular  potential 

well  depth, (C/k)  K 

nc 

noncondensable  molecxilar  weight.W^^  Ibm/lbm-mole 
condensable  molecular  weight,  Ibm/lbm-mole 

total  static  pressure,  p  psia 

partial  pressure  of  vapor,  p^  psia 

Output  -  gas  mixture  viscoMty,  /x  Ibm/hr-ft 

Limits  -  individual  gas  mixture  component  viscosities  are  given  by 
equation  8.2-18  of  Reference  23.  Gas  mixture  viscosity  is  given  by 
equation  8.2-30  of  Reference  23.  The  binary  diffusion  coefficient  is  given 
by  equation  8.2-44  of  Reference  23.  Molecular  collision  diameter  and 
molecular  potential  well  depth  of  the  mixture  are  given  by  equations  8.  6-3 
8.  6-4  and  8.  6-5  of  Reference  23.  Collision  integrals  are  approximated 
from  Table  I-M  of  Reference  23. 

Reference  -  Reference  23. 

GASP4: 

Description  -  Calculation  to  determine  viscosity  of  vapor  at  a  given 


temperature. 
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Input  -  temperature,  T  K 

condensable  molecular  collision 

diameter.  O  A 

V 

condensable  molecular. potential 

well  depth,  (C/k)^  K 

condensable  molecular  weight,  Ibm/lbm-mole 

Output  -  viscosity,  Ibm/hr-ft 

Limits  -  viscosity  is  given  by  equation  8.2-18  of  Reference  23.  Collision 

integrals  are  approximated  from  Table  I-M  of  Referehce"23. 

Reference  -  Reference  23. 


FORCE; 

Description  -  calctdation  to  determine  phase  interface  resultant  friction 
force  and  mixing  section  wall  resultant  friction  force  for  a  giyen  mixing  sec¬ 
tion  and  liquid  jet  geometry  and  known  Reynolds  numbers. 

Input  -  average  liquid  jet  diameter,,  D,  in. 

L 

average  mixing  sectipn  diameter,  D  jin. 

8 

niixing  section  length,  L^^ 

+  average  bulk  liquid  velocity,  V..  ft/sec 

average  bulk  gas  velocity,  V  ft/sec 

8 

+  average  relative  velocity,  V  ft/sec 

8  2 

liquid  jet  surface  area,  (S.)22 

bulk  gas  Reynolds  number  at 
state  1,  (Re)  , 

81 

bulk  gas  Reynolds  number  at 

state  2,  (Re)  _ 
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+  For  mixing  section  Analytical  Models  I,  Ila.  and  III  the  phase  interface 
friction  shear  stress  is  based  on  the  relative  velocity,  V'  ,  and  for  Analytical 

Model  lib  it  is  based  on  the  velocity  (Vg-Vj^). 
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Reynolds  nvmber  of  anniilar  gas  stream 

at  state  (Re)  . 

rgl 

Reynolds  nximber  of  annular  gas  stream 


at  state  2,  (Re)  . 

rg2 

3 

gas  density  at  state  1, 

Ibm/ft 

3 

gas  density  at  state  2,  p  ^ 

Ibm/ft 

-  phase  interface  resultant  friction 

force ,  F. 

1 

-  Ibf 

mixing  section  wall  resultant  fric- 

tion  force,  F 

Ibf 

w 

References  -  friction  factors,  f.  and  f  ,  from  equations  E32-E35. 

:  ^  1  w  2 

Shear  stress  from  equation  E31.  Shear  forces  from  equations  E29  and 
E30. 
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Static  Pressure  at  Mixing  Section  Inlet,  Pj(psia) 

FIGURE  3  -  LIQUID  JET  BREAKUP  LENGTH  VERSUS 
STATIC  PRESSURE  AT  MIXING  SECTION  INLET  - 
CONDENSABLE  CASE 
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FIGUR£  4  •  Ratio  of  Final  Mean  Drop  I 
Radius  to  Initial  Jet  Radius  versus  Mask 
\  Flour  Rate  Ratio  •  Noncondensable  Clf  e 
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FIGURE  5  •  Ratio  of  Final  Msan  Drop 
Radius  to  Initial  Jet  Radius  versus  Mass 
Flow  Rate  Ratio  •  Condensable  Case 


10.  5. 


1.  0.5 


Mass  Flow  Rate  Ratio,  w  «  (w.  /w  ) 

^  go' 


(eqn.  A 10) 
(eqn.  All) 

(eqn.  A9) 


Ratio  of  Final  Mean  Drop  Radio*  to  Initial  Jet  Radius.  — 


r^y  calculated  from  eqn.  A9 


*ol 

I 


Pj.P*»a 


Test  results  of  \ 
RUN  NO.  1 
(honcondensable  flow) 


■  323.  Ibm/min 
w  •  13,0  Ibm/min 

VO 

w  ■  2.  52  Ibm/min 

nc 

T  a  340.  F 

go 

T_  a  40.  F 
Lo 

Pgo  ■94.5p#ia 

Plo  • 

fa>  a  20.81 

o 

■  0,  147  in. 

t-  a  0.  651  in 


Pj.psia 


Test  results  oOf 
RUN  NO.  2  ^^0 

(condensable  flow)  \ 


0.05 


a  358.  Ibm 
w  a  23.  5  Ibm 

VO 

T  a  283.  F  “ 

go 

T_  a  40.  F 
Lo 

P„_  ■  36.  psia 
go 

Plo  * 

w  a  15.23  _ 

o 

^L1  * 

■  0. 651  in. 


0.01 


0.005 


FIGURE  6  -  Volumetric  Flow  Ra^  Ratio  of  Liquid 
Phase  to  Gas  Phai&e  versus  Ratio  of  Final  Mean  Drop 
Radius  to  Initial  Jet  Radius  for  a  Condensable  and 
Noncondensable  Case 


0.001 


1,  0.5  O.'i  0.05  0.01  0.  005  0. 

Volumetric  Flow  Rate  Ratio  of  Liquid  Phase  to  Gas  Phase,  w 

vol 


Mixing  Section  Temperature*,  (F) 


101 


FIGURE  7  -  VARIOUS  MIXING  SECTION  TEMPERATURES 
VERSUS  DISTANCE  FROM  MIXING  SECTION  INLET 
PLANE  -  MODEL  I 
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FIGURE  8  -  VARIOUS  MIXING  SECTION  VELOCITIES 
VERSUS  DISTANCE  FROM  MIXING  SECTION  INLET 
PLANE  -  MODEL  I 


Transfer  Coefficients.  Btu/hr-ft 


p  ■  36.  0  psia 

go 

PLo  “^^^.Opeia 
p  s30.  Opeia 


0. 295  in 
283.  0  F 


T_  e40.0F 
Lo 


15.44 


No.  of  Jets  •  1 


Constant  pressure  mixing,  no  none  onde  ns  able 


500, 00( 


bulk  gas  mixture  Reynolds  no.  •  (Re) 


bulk  gas  mixture  Reynolds  no.  relative  to 
liquid  jet  interlace,  (Re) 


'  v 

bulk  liquid  jet  Reynolds  na  relative  to  liquid  .  \ 
jet  interface,  ' (Re)  'v 


heat  transfer  coefficient  from  liquid  jet  interface 
to  bulk  liquid,  h 


overall  heat  transfer  coefficient,  U 


im 


heat  transfer  coefficient  from  bulk  gas  mixture 
to  liquid  jet  interface,  h-  I 


Distance  from  Mixing  Section  Inlet  Plane,  L{in) 


FIGURE  9  -  VARIOUS  REYNOLDS  NUMBERS  AND  HEAT 
TRANSFER  COEFFICIENTS  VERSUS  DISTANCE  FROM 
MIXING  SECTION  INLET  PLANE  -  MODEL  I 
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FIGURE  U  -  VARIOUS  MIXING  SECTION  TEMPERATURES 
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FIGURE  14  -  LIQUID  JET  AND  MIXING  SECTION 
DIAMETERS  VERSUS  DISTANCE  FROM  MIXING 
SECTION  INLET  PLANE  -  MODEL  Ila 
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FIGURE  15  -  VARIOUS  MIXING  SECTION  TEMPERATURES 
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SECTION  INLET  PLANE  -  MODEL  Ilb 
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FIGURE  19  -  VARIOUS  MIXING  SECTION  TEMPERATURES 
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MODEL  ni 


Mixii^  section  V«locitiM ,  (ft/tecj' 


il4 


V 

*  ^IV5  pui* 

T 

id 

»  340.  0  F 

■  2L25 

? 

■  pfia  ' 

^X,6 

•  40/6  F 

6  ■  0.  338 

^  p 

Pi 

•  «?..  7  ppia 

•^1,1 

•  0.  245  in. 

Np,  of  JpUi  ■  1.0 
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PLANE  -  MODEL  lU 
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FIGURE  21  -  VARIOUS  REYNOLDS  NUMBERS,  HEAT  TRANSFER 
COEFFICIENTS,  AND  MASS  TRANSFER  COEFFICIENT  VERSUS 
DISTANCE  FROM  MIXING  SECTION  INLET  PLANE  -  MODEL  HI 


Raynolds  Numbara 


Liquid  Jet  and  Mixing  Section  Diameters ,  (in) 


116 


p _ 

«  94. 5  psia 

T 

*  340.  0  F 

o>  =21.25 

•« 

go 

*  40. 0  F 

o 

*  325.0  psia 

6  =0. 838 

^  o 

•*1 

■  89. 7  psia 

=  0.295  in 

No.  of  Jets  =  1.0 

FIGURE  22  -  LIQUID  JET  DIAMETER,  MIXING  SECTION 
DIAMETER,  AND  TOTAL  STATIC  PRESSURE  VERSUS 
DISTANCE  FROM  MIXING  SECTION  INLET  PLANE  - 
MODEL  m 


Mixing  Section  Total  Static  Pressure,  p  (psia) 


P«rc«nt,ol  Vapor  Condenscid 


117 


Plo 

-  325.  0  psia 

^Lc 

=  40.  0  F 

4  =  0,  8S8 

^  o 

Pi 

°L1 

*  0.  295  in 

No.  of  Jots  s  1.  0 

T 

go 

s  500.  0  F 

u 

o 

*  21.24 

FIGURE  23  -  PERCENT  OF  VAPOR  CONDENSED  VERSUS 
DISTANCE  FROM  MIXING  SECTIO.^  INLET  PLANE  FOR 
VARIOUS  INLET  MACH  NUMBERS  -  MODEL  Ila 


Per  cant  of  Vapor  Condensed 


■  94. ‘5  psia 

T 

*  340.  0  F 

go 

«  325.  0  psla 

*  40.  0  F 

■  90.  0  psia 

®L1 

•  0.295  F 

o>  *21.24 
o 

No.  of  Jets  *1.0 


Constant  pressure  mixing*  with  noncondensable 


Mixing  Section  Diameter,  D  (in) 


120 


Input  Data; 

Model  Reference 


I  Figure*  7>10 

He  Flguree  11«14  No.  of  Jets  «  1 . 0 

nb  Figures  15-18 


Distance  from  Mixing  Section  Inlet  Plane,  L(in) 
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